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Abstract 



Gravitational wave signatures from cosmic strings are analyzed numerically. 
Cosmic string networks form during phase transistions in the early universe 
and these networks of long cosmic strings break into loops that radiate energy 
in the form of gravitational waves until they decay. The gravitational waves 
come in the form of harmonic modes from individual string loops, a "confu- 
sion noise" from galactic loops, and a stochastic background of gravitational 
waves from a network of loops. In this study string loops of larger size a 
and lower string tensions (where is the mass per unit length of the 
string) are investigated than in previous studies. Several detectors are cur- 
rently searching for gravitational waves and a space based satellite, the Laser 
Interferometer Space Antenna (LISA), is in the final stages of pre- flight. The 
results for large loop sizes (a = 0.1) put an upper limit of about G^i < 10~^ 
and indicate that gravitational waves from string loops down to G// ~ lO"^*^ 
could be detcctabc by LISA. The string tension is related to the energy scale 
of the phase transition and the Planck mass via G/x = A^/m^j, so the limits 
on Gfi set the energy scale of any phase transition to < lQ~^-^mpi. Our 
results indicate that loops may form a significant gravitational wave signal, 
even for string tensions too low to have larger cosmological effects. 
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small loops |86 



7 



7.4 The gravitational wave energy density per log frequency from 
cosmic strings for a = 10~^ and 7 = 50. Smaller a reduces 
the background for equivalent string tensions, with a LISA 
sensitivity limit of G/x > 10~^^ if the Sagnac sensitivity limits 
are used. The bars below the main LISA sensitivity curve are 
the Sagnac improvements to the sensitivity, while the dotted 
line is the Galactic white dwarf binary background (GCWD- 
WD) from A.J. Farmer and E.S. Phinney, and G. Nelemans, 
et al., [371 [76]. The dashed lines are the optimistic (top) and 
pessimistic (bottom) plots for the extra-Galactic white dwarf 
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Chapter 1 
Introduction 



1.1 Cosmic Strings 

Quantum field theory predicts that the universe can transition between dif- 
ferent vacuum states while cooling during expansion. In these theories the 
vacuum expectation value of a field (0) can take on non-zero values when 
the ground state, at the minimum of the zero-temperature potential or "true 
vacuum" , breaks the symmetry of the underlying Lagrangian [1061 HI] . In 
the hot early universe, thermal effects lead to an effective potential with a 
temperature dependence, such that the vacuum starts in a symmetric "false" 
vacuum and only later cools to its ground state below a certain critical tem- 
perature [66l 11061 [HI [981 1104] . In general however the expansion of the uni- 
verse occurs too quickly for the system to find its true ground state at all 
points in space [57]. 

Depending upon the topology of the manifold of degenerate vacua, topo- 
logically stable defects form in this process, such as domain walls, cosmic 
strings, and monopoles [1061 HH [9H1 I104[ [85] . In particular, any field the- 
ory with a broken U(l) symmetry will have classical solutions extended in 
one dimension, and in cosmology these structures generically form a cos- 
mic network of macroscopic, quasi-stable strings that steadily unravels but 
survives to the present day, losing energy primarily by gravitational ra- 
diation [103[ HH [HHl [m inn [S]- A dual superstring description for this 
physics is given in terms of one- dimensional branes, such as D-strings and 
F-strings [1051 EQl |35l ESI [Ml E] ■ 

Previous work has shown that strings form on scales larger than the 
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horizon, and expand with the horizon with cosmic time [HS]. These are 
referred to as "infinite strings" , and through various mechanisms give rise to 
loops which form as a fraction of the horizon. These loops are the primary 
potential source of gravitational radiation detectable from the earth. 

Calculations are undertaken at much lower string masses than previously, 
with several motivations: 

1. Recent advances in millisecond pulsar timing have reached new levels 
of precision and are providing better limits on low frequency back- 
grounds; the calculations presented here provide a precise connection 
between the background limits and fundamental theories of strings and 
inflation [391 [53l [551 [97] . 

2. The Laser Interferometer Space Antenna (LISA) will provide much 
more sensitive limits over a broad band around millihertz frequencies. 
Calculations have not previously been made for theories in this band of 
sensitivity and frequency, and are needed since the background spec- 
trum depends significantly on string mass [TOl [H [22] . 

3. Recent studies of string network behavior strongly suggest (though they 
have not yet proven definitively) a high rate of formation of stable 
string loops comparable to the size of the cosmic horizon. This results 
in a higher net production of gravitational waves since the loops of 
a given size, forming earlier, have a higher space density. The more 
intense background means experiments are sensitive to lower string 
masses [gilfTIMITT]. 

4. Extending the observational probes to the light string regime is an 
important constraint on field theories and superstring cosmology far 
below the Planck scale. The current calculation provides a quantitative 
bridge between the parameters of the fundamental theory (especially, 
the string tension) , and the properties of the observable background [581 
[85]. 

Cosmic strings' astrophysical properties are strongly dependent upon two 
parameters: the dimensionless string tension G/i (in Planck units where c = 1 
and G = mpC"^) or G^ijc? in SI units, and the interchange probability p. 
Our main conclusion is that the current pulsar data [531 iSSl [90] already 
place far tighter constraints on string tension than other arguments, such as 
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microwave background anisotropy or gravitational lensing. The millsecond 
pulsar limits are dependant upon the nature of the source, and are detailed 
in later chapters. Recent observations from WMAP and SDSS [109] have put 
the value of G/i < 3.5 x 10~^ due to lack of cosmic background anisotropy 
or structure formation consistent with heavier cosmic strings. From [109] it 
is found that up to a maximum of 7% (at the 68% confidence level) of the 
microwave anisotropy can be cosmic strings. In other words, for strings light 
enough to be consistent with current pulsar limits, there is no observable 
effect other than their gravitational waves. The limit is already a powerful 
constraint on superstring and field theory cosmologies. In the future, LISA 
will improve this limit by many orders of magnitude. 

The energy density of the cosmic strings is dependant upon the energy 
scale of the associated phase transition: for a GUT (Grand Unified Theory) 
scale phase transistion the energy is of the order 10^^ GeV and the elec- 
troweak scale on the order of 10^ GeV. In the case of cosmic strings, this 
energy determines the mass (energy) per unit length fi of the string, and the 
dimensionless string tensions Gfi/c^ in SI units or G/j, when c = 1. For this 
study lighter strings on the order of 10~^^ < Gfi < 10~^ are used, as analysis 
of millisecond pulsars has ruled out higher mass strings [53] . 

In general, lower limits on the string tension constrain theories farther 
below the Planck scale. In field theories predicting string formation the 
string tension is related to the energy scale of the theory through the 
relation Gfi oc A^/m^; [1061 IFTl [H]. Our limits put an upper limit on A^ in 
Planck masses given by A^ < lO""^'^, already in the regime associated with 
Grand Unification; future sensitivity from LISA will reach A^ ^ 10~^, a 
range often associated with a Peccei-Quinn scale, inflationary reheating, or 
supersymmetric B-L breaking scales [52] . In the dual superstring view, some 
current brane cosmologies predict that the string tension will lie in the range 
10~^ < Gfi < 10~^^ [961 [Ml [97]; our limits are already in the predicted range 
and LISA'S sensitivity will reach beyond their lower bound. 

1.2 Equation of State and Thermodynamics 

In this section a selection of the constituents of the universe and some of 
their thermal properties are analyzed P[T9l[33l[lIl[65l[73l[71l[IHl [791 [T07]. 
These results are important in calculations of the rate of cosmic expansion 
which, in turn, affects the number of cosmic string loops formed and their 
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size. For a more thorough discussion of relativistic cosmology and general 
relativity see Appendices A and B and the aforementioned references. 

We can define a region of space as a comoving volume denoted by a{t)^ 
where a{t) is the scale factor, which expands with the universe. Of interest 
in this region are the equations of state of the constituents and how they 
evolve with time. Of utmost importance to cosmic dynamics is the stress- 
energy tensor at different points in space. More details on the stress-energy 
tensor are given in Appendices A and B, but it is a second rank tensor that 
describes the mass/energy content, energy flux, momentum density, and stess 
at a spacetime point. 

As a useful and relevant example, a perfect fluid has a stress-energy tensor 
of the form, 

T^,u = {p + p)Uf,Uu + pg^u, (1-1) 

where u is the four velocity of the fluid. In a local (Minkowski) frame the 
tensor is: 

Too = p and Tij = p6ij, (1.2) 

where p is the density and p is the pressure of the cosmic constituent. In this 
frame there is a simple definition of the stress-energy tensor: the 00 com- 
ponent is the energy density and the diagonal components are the pressure. 
Again, see Appendices A and B for more details. 

In general we can relate the equation of state to the energy density p via, 

p — wp. (1.3) 

Conservation of the stress-energy tensor gives a relationship between the scale 
factor and density, 

d^T"" + T';;^T^'' + ri^T^'^ = 0, (1.5) 

where F^^ are the Christoffel symbols (see Appendix A) and is the covari- 
ant derivate. Using the p — component of the stress-energy tensor leads 
to, 

p + 3(p + p)- = 0. (1.6) 
a 

This equation gives the time dependance of the density of each constituent, 
and also describes the dependance of the densities on the scale factor (or 
volume) . 
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For the constituents of the universe we have: 



Non-relativistic matter, w = 0, p (X a 
Ultra-relativistic matter, w = -, p oc a 



-3 



(1.7) 
(1.8) 



-4 



Vacuum energy (Dark energy), w = —1, p = constant. (1-9) 

These relations indicate a universe composed of "dark enegy" eventually has 
the expansion dominated by it. Both relativistic and non-relativistic matter 
density shrink as the universe expands, which corresponds to the observation 
that the very early universe was dominated by radiation, which was replaced 
by non-relativistic matter, finally to be overtaken today by vacuum energy. 
A more complete description is given in Appendix B. 

For reference we investigate the long "infinte" string equation of state |18] , 



which varies with the rms speed of the long strings. Assuming the infinite 
strings are ultra-relativistic we find. 



For the dynamics of the scale factor, the cosmic strings are a miniscule 
component, thus their equation of state is unimportant. The upper limit 
on the density of strings is very small from observations of the cosmic mi- 
crowave background looking for string related inhomogeneities using WMAP 



and SDSS ^109j . 

1.3 Field Theory and Symmetry Breaking 



The Lagrangian density C defines the spectrum of particles and interactions 
in a theory. A number of requirements are imposed on any viable field the- 
ory to include: the theory must be Lorentz invariant (relativistic), invariant 
under general gauge transformations, and result in a spectrum of particles 
and interactions that matches observation. Said theory must also take into 
account local symmetries associated with the gauge group of the field. For 
example, the electromagnetic interaction is symmetric under the U{l)em sym- 
metry group, while the electroweak interaction is symmetric under the gauge 
group SU{2)l X f/(l)y. QCD is symmetric under SU{?>)qcd- 





_ 1 
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For example, at some finite temperature the symmetry of tlie SU{2)l x 
U{1)y group is broken to leave the f/(l)e,„ unbroken. The resulting spectrum 
is the one in which we are familiar: spin-1 photons and spin-1/2 leptons. This 
is predicted to occur at a temperature of about 10^ GeV, or T ^ lO^^K, using 
k = 8.617 X 10-^ GeV/K. Likewise, the group SU{3)qcd x SU{2)l x U{1)y 
can be an unbroken group of another phase transition, e.g. the SU{5) group. 

For more detail let us look at the electroweak model. The Langrangian 
for the electroweak theory has a pure gauge and a Higgs boson part, given 
by: 



jO^ew = fil^D^ - m)f + Cgauge + ^Higgs, (1-12) 



where 



Cmggs = {D,4>)\D^4>)-V{4>\4>.T). (1.14) 



Cgauge = —W^^.W^ - -B ,,B>^^ , (1.13) 



and 



= d,W: - d.W; - ge'^'^'Wlw;, (1.15) 
B^. = d^B,-d,B^ (1.16) 

The are the three gauge fields associated with the generators of SU (2) 
and B^j, is the guage field associated with C/(l)y: Here is the covariant 
derivative and is 

D,^d, + ^Tj¥;+'-^B, (1.17) 

where the are the three generators of the group, represented by Pauli 
matricies for SU{2). g and g' are the gauge couphng constants for SU{2) 
and f/(l). 

The potential V{(f),T) determines the behavior of the theory at varying 
temperatures, T. At high T the potential keeps the symmetry of the under- 
lying theory. When T drops to a critical value the old symmetry is broken, 
and the spectrum of particles is changed. For the electroweak theory, the 
massless exchange particles, the W and B, which act like a photon, are gone 
and replace by the Z and gauge bosons, which are massive. The electro- 
magnetic field is massless, from the superposition given by, 

A^^cos{9w)B^ + sm{9w)W^, (1.18) 



16 



which remains massless after the break to f/(l)em symmetry. The 9w is the 
Weinberg angle, which when coupled to the electromagnetic interaction with 
strength e gives, 

e = g sm{9w) = g' cos{9w) ■ (1.19) 

An elementary example of this symmetry breaking is the "mexican hat" 
potential. At the critical temperature the minimum of the potential V"(</)) 
takes on a set of degenerate values that break the symmetry of the theory. 

Theories with global symmetry will tend to radiate Goldstone bosons, 
which typically reduce the lifetime of string loops to approximately 20 pe- 
riods. This reduces to possible gravitational wave signal, likely making it 
invisible. Strings created by a locally symmetric gauge theory that radiate 
only gravitational waves are analyzed here, as most of the energy lost from 
the strings is in the form of gravitational waves. 



1.4 Vacuum Expectation Value of Fields and 
Vacuum Energy 

Of great importance when dealing with topological defects is the vacuum 
expectation value (VEV) of the field. The VEV is the average value one 
expects of the field in the vacuum, denoted by (0). This can akin to the 
potential: if the potential is 100 V throughout space that has no physical 
meaning. But in the case of a homogeneous real scalar field, the VEV can 
lead to a non zero vacuum energy given by the stress-energy tensor, 

dC 

which can be solved for the 00 term to get, 

To^^p=\[^ + {V<t>f]+V{(t>). (1.21) 

From above we find that even if (0) is stationary, homogenous, and zero 
the potential V{{(j))) can result in a vacuum energy density. 

During a phase transition the old VEV of zero changes to some non-zero 
value. Given the finite travel time of information, for first order phase tran- 
sitions bubbles of this new vacuum propagate at the speed of light. Where 
these bubbles intersect can lead to cosmic turbulence as well as topological 
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defects, to include cosmic strings. For our strings, the zero VEV region is 
within the string and the new VEV without. 

The regions of old vacuum (topological defects) are stable if the thermal 
fluctuations in the field are less than the differences in energy between re- 
gions, i.e. the old and new vacuums. This is typically given by the Ginzburg 
temperature Tq and is close to the critical temperature of the phase tran- 
sition. Below the Ginzburg temperature the defects "freeze-out" and are 
stable [T06] . 



1.5 Cosmic String Motion 

Cosmic string motion can be characterized by paramaterization of the four 
space-time dimensional motion onto a two-dimensional surface called a world- 
sheet. With the coordinate parameterization given by: 

x^" = X^'ia''), a = 0,1, (1.22) 

where the x'^ are the four spacetime coordinates, X^^ are the function of a", 
and the a" are the coordinates on the worldsheet of the string. In general 0"° 
is considered timelike and is considered spacelike. On the worldsheet the 
spacetime distance is given by, 

ds^ = g^^dx^dx", (1.23) 

= ^abda''da\ (1.25) 

where 7afe is the induced metric on the worldsheet. 

The equations of motion for the string loops devised by Yoichiro Nambu 
were formulated by him assuming that that the action was a functional of 
the metric g^i, and on functions X^. Also added were the assumptions: 

• the action is invariant with respect to spacetime and gauge transfor- 
mations, 

• the action contains only first order and lower derivatives of X^. 
Ths solution is dubbed the the Nambu-Goto action: 

S = -fi ^y^d^a, (1.26) 
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where 7 is the determinant of '~)ab- This is not the only relativistic action for 
cosmic strings, but is the only one that obeys the caveats above. Including 
higher order derivatives corrects for the curvature of the string, but is only 
important in the neighborhood of cusps and kinks ^ 
Minimizing the action leads to the following equations, 

^ +r:,— ^ = (1.27) 

where F are the Christoffel symbols. 

In flat spacetime with g^^ = and F^^ = the string equations become, 

^ ^ ^r'^) = 0. (1.28) 



da"' \ da^ 

Now we use the fact that the Nambu action is invariant under general 
coordinate reparameterizations, and we choose a specific gauge. This does 
not disturb the generality of the solutions we find, since the action is invariant 
under these gauge transformations. Here we choose: 

701 = 0, (1.29) 
700 + 711 = 0, (1.30) 

this called the "conformal" gauge, because the worldsheet metric becomes 
conformally flat. 

Our final choice is to set the coordinates by 

t = x° = a° (1.31) 
a = a\ (1.32) 

so that we have fixed the time direction to match both our four dimen- 
sional spacetime and the worldsheet. This is sometimes called the "temporal" 
gauge. This leads to the equations of motion, 

i-x' = 0, (1.33) 
x2 + x'2 = 1, (1.34) 
x-x" = 0. (1.35) 

The first equation states the orthogonality of the motion in these coordi- 
nates, which leads to x being the physical velocity of the string. The second 
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equation is statment of causality as the total speed of the motion is c, which 
can also be interpreted as, 



da = —j^=^= = dE/ii, (1.36) 

So the coordinates a measure the energy of the string. 

The final equation is a three dimensional wave equation. This equation 
determines the behavior of the strings and has a general solution of the form: 

x(cx,t) = i[a((7-t) + b((7-t)], (1.38) 

where constraints from the equations of motion give, 

a' = b' = 1. (1.39) 

We can also write down the stress-energy tensor, momentum, and angular 
momentum in this choice of gauge: 

r-^,/<i.[i^i"-.V|^«(x-x(M)). (1.40) 

which leads tot he energy given above, as well as, 

p = H J da±{a,t), (1.41) 

J = II J dax{a,t) x±{a,t). (1.42) 

A cosmic string network is formed and behaves as a random walk. As 
the long strings and string segments move and oscillate they often intersect 
one another. Depending upon the characteristics of the strings a loop can be 
formed, which we discuss in the next section. 



1.6 Cosmic String Loops 

Previous results have shown that the majority of energy radiated by a cos- 
mic string network comes from loops which have detached from the infinite 
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Figure 1.1: Loop formation from long strings. Two possibilities for the cre- 
ation of cosmic string loops from the intersection of long strings. 

strings. Loops are formed by interactions between the infinite strings, as 
shown in Fig. II. 1[ Of utmost importance is the interchange probability of 
the strings p. If p=l then each time strings intersect they exchange partners, 
and a loop is much more likely to form. 

Numerical simulations, analytic studies, and heuristic arguments indicate 
cosmic string loops form on the "one-scale" model; that is, their size scales 
with the horizon as in Fig. 11.21 The ratio of the size of the loops to the 
horizon is given by a, and ranges from the largest size of 0.1 down to the 
miniscule 10~^. There is some debate about the stability of the larger string 
loops, and older numerical simulations predict that any large loops should 
decay to the smaller sizes. More recent simulations indicate the possibility of 
large stable loops, which is part of the motivation for this study I1U2[ [7T] . 

To analyze the strings in the universe we must make estimates of the 
loop formation rate, based on our models of string creation. First we analyze 
the number of infinite strings created within the cosmic horizon. By infinite 
we mean they extend beyond the horizon, and stretch as the horizon grows. 
From this the number of loops formed per unit time can be analyzed. 

To do this we define three parameters: A is the number of infinite strings 
within the horizon, B is the creation rate of loops within the horizon, and 
a is the size of loops formed relative to the horizon. For the first two we 
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Figure 1.2: Scaling of cosmic string loops. The figure on the left shows 
a newly formed loop, scaled to a certain fraction of the horizon. When 
the horizon expands as on the right, the newly formed loop on the right 
scales with the expansion, while the old loop on the left continues to radiate 
gravitational waves and shrink. 



have [18], 

A = p^^, (1.43) 

(1.45) 

where poo is the energy density of infinite strings. In principle the num- 
ber density at any time can be calculated after these parameters have been 
chosen. 

For this study we set the number density of loops per log cosmic time 
interval in years of 0.1 , 



a \ c J \a[to) J 

where tc is the time of creation of the loops and to is the time of observation. 
In this equation our paramter A^^ has replaced A and B and incorporated 
Poo- Results are discussed in more detail in the chapters describing the grav- 
itational wave signal results. 
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Chapter 2 

Introduction to Gravitational 
Waves 



Einstein's theory of general relativity is the foundation of modern cosmology, 
but at its inception it was a theory that seemed to be un-needed. Today, 
though, it is found indespensable in explaining a large number of phenom- 
ena and in the proper functioning of such things as the Global Positioning 
System. One of the hallmarks of general relativity, and indeed of all metric 
theories of gravity, is the production of gravitational waves. The motion of 
a gravitational mass alters the space-time nearby, and this alteration then 
propagates at the speed of light producing a gravitational wave. 

In this chapter we discuss the basics of gravitational waves and introduce 
some of the equations that are used to describe them. Of great import 
are the plane wave solutions, given that most of the sources are a distance 
from the detector. For more details on general relativity see Appendix A 
and references cited there. For more on gravitational waves in particular 
see [221 [70]. 



Einstein's field equations are given by the relationships between the curvature 
of spacetime to the stress-energy tensor: 



2.1 Einstein' 



s Equations 




(2.1) 



where Rf^u is the Ricci tensor and R is the Ricci scalar g^^R'^'-' 
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For the weak field case we assume a small perturbation on our Minkowski 
metric so 

g^lu = Vf^u + (2.2) 
where h^y << 1. The linearized Ricci tensor then becomes 

= daV^^ + duT'^^, (2.3) 

where the F^^ are Christoffel symbols. 

After inserting this into Einstein's equations we get: 

- V - V^ud'^d^Kf^ + d^d'^K^ = levrGT^,, (2.4) 

where h^^ = h^y — \rj^yh is the trace- reversed tensor. 

2.2 Lorentz Gauge and Wave Equation 

From the above equations it isn't clear there are wave solutions, but if we 
choose the gauge 

5" V = 0, (2.5) 
which is known as the Lorentz gauge, we find the equation is now, 

= -levrGT^,, (2.6) 
(-92 + V2)V = -167rGr^,. (2.7) 

If we take the trace of both sides of the equation we find dad°'h = IGvrG T. 
If this is inserted into the above equation we get 

i-d^ + V^) V = -ISttG (t,, - ^r]^,T^ , (2.8) 

where T = g^^T^". Note that Eqs. 12.71 and 1 2.81 are equivalent, although in 
practice Eq. 12.81 is more useful when the stress-energy tensor is known. This 
looks like a traditional wave equation, with the only difference being the 
tensor nature of the spacetime metric. The general solution to Eq. 12.81 is 
given by the retarded time solution, 

V(x, t)=AG ! ^^^ M^^ t-\^ -^\) ^ ^2.9) 
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Of course now we look at a region of space devoid of any matter, in which 
case the Einstein equations become, 

(-a2 + v2)V = o. (2.10) 

We should expect to get plane wave solutions from this. 



2.3 Wave Solution and TT Gauge 



2.3.1 Plane Waves 

So solutions will be of the form 



h 



(2.11) 



where A^^iy is a tensor describing the polarization of the plane waves. 

The Transervse Traceless or TT gauge is used to greatly simplify the 
appearance of the wave solutions to Einstein's field equations. The generality 
of the results is not damaged by this particular gauge choice, so it is to our 
benefit to make this change. 

There are three main steps to creating the TT gauge: (a) from the Lorentz 
gauge we find that A^^k^ = and /c^/c^ = 0. (b) We can also choose a 
Lorentz frame in which our velocity u is orthogonal to A so that Af^yu'^ = 0. 
(c) Finally we choose tr{A) — A'^^ — 0. 

If we pick our particular frame to have the four velocity given by Ufj, — 
(1,0,0,0), 





h 



/o 








K 
K 












0/ 



-i{kz—wt) 



This can be rewritten in terms of polarization tensors e 



IXV-i 



a e'fe-'^''^-'"^^ + be 



P-f -i{kz—wt) 
X ^ 



(2.12) 



where a and b are complex constants. 



25 



2.3.2 Spherical Waves 

For spherical waves we can use the same ansatz as before, but with 1/r 
dependence added: 

= (2.13) 

r 

In this case we make sure we have the proper sign in our exponential to 
ensure an outgoing wave. 



2.4 Particle Orbits from Gravitational Waves 

Now the important question: what happens to a particle as the wave passes? 
Nothing! Let's clarify that statement somewhat, lest the reader assume the 
author suggests there is no visible effect of a passing gravitational wave. From 
the geodesic equation 

— = -TlXu\ (2.14) 

using the TT gauge we find du'^/ dr = and the particle remains at rest in its 
reference frame. So in fact there is no observable effect on a single particle. 
This is due to our choice of gauge, so a different approach must be taken. 
We must look at two or more particles to observe an effect. 

More specifically, we want to look at the variation of the geodesies of two 
particles as a gravitational wave passes. For two geodesies, one at x and the 
other ai X + X their deviations can be measured by 



^ = (2.15) 

where X" is the four vector separation between the geodesies, and we use 
the definition: 

dX"" dr" 
DX^/Dt = X^^^^ + (2-16) 

Only the R^q^q terms will contibute due to our TT gauge choice. 

In reality we must look not at the TT gauge, but at an actuall laboratory 
frame. This is denoted the proper lab frame and to first order in h we find 
that T — t, so we can write, 

^ = -c'RVoX\ (2.17) 
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Figure 2.1: Diagrams showing the motion of a ring of particles under the 
influence of a passing gravitational wave. The top is plus polarization and 
the bottom is cross polarization. 



which results in the differential equation, 



This can be solved to give us 

AX^ = hj^XK (2.19) 

Here the X^ are the three vector components of the coordinate separation 
vector. We see that the oscillations of spacetime are orthogonal to the direc- 
tion of propogration of the wave. 

In greater detail, the motion of the particles is determined by the po- 
larization of the wave: "cross" and "plus" polarizations produce different 
relative motions. Fig. 12.11 shows the motion of a group of test particles when 
a gravitational waves is incident upon them. 
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2.5 Characteristic Strain 

In the transverse traceless gauge we can define the wave in terms of its Fourier 
components by, 

/oo p /'27r 

df / dk / d^i^(K{fXi^)etj{ki^) + 
oo .In Jo 

+ K if, k, ij)e^j{k, V^))e-^2'^^(*-^-^), (2.20) 

where, hA{—f,k,ip) = h\{f^k^ip). Also ^ is a unit vector pointed in the 
direction of propagation and ■0 is the angle of polarization of the wave. The 
two polarization tensors have the orthogonality relation e4(A;)e^,(fc) = 25'^^ , 
no sum over A. Now we define the spectral density Sh{f) from the ensemble 
average of the fourier amplitudes, 

h:,{f, k, ^)hA'(f', k', V'O) = SAA'Sif - f) 

where 6^{k, k') = 5(0 - 0')5(cos(^) - cos(^'))- 

Now we write down the average of the real waves and the result is. 



oo 



(/i,,(t,x)/i*^(t,x)> = 2 / dfSuif), 

J — oo 

= 4/ di\nf)fSM, (2.22) 

Jf=0 



where we have used a one-sided function for 3^.. 
The characteristic strain is defined by. 



{hij{t,^)h'^{t,^)) = 2 / ci(ln/) hlif). (2.23) 

Jf=0 

Therefore, 



hl{f)^2fSM, (2.24) 

(2.25) 
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which is an important result used in the analysis of spectra from a gravita- 
tional wave source. 

In the normalization used in the LISA sensitivity plots [63j we get the 
relation, 

SH{f) = \Hf)f . (2.26) 

There are numerous definitions of strain in the multitude of publications 
on the subject of gravitational waves. The author wishes to maintain the 
notation as close as possible to the literature, and to be consistent throughout 
this work. To that end, the above definitions are adhered to in all that follows. 

We have derived the basic formulas for gravitational waves from linearized 
general relativity. Of note is the use of the transervse, traceless (TT) gauge, 
in which the waves are transerve and are superpositions of two polarizations, 
+ and X . A single particle feels no acceleration from the passing wave as its 
coordinate system moves with the waves, but several particles will observe 
relative accelerations between them. 



29 



Chapter 3 



Power Radiated in 
Gravitational Waves 



The universe is rich in sources of gravitational waves; from black hole binaries 
to cosmological phase transitions, the spectrum should be quite full. The 
keys to detection are the power radiated by the source, its distance from 
the observer, and the nature of the signal. In this chapter we overview the 
primary equations of the power radiated by a source, and discuss a few of 
the sources. For a more detailed survey readers should consult [221 HOI 191] . 

The linearized solutions to Einstein's equations are expanded in the re- 
tarded time for a low velocity source, and terms are collected into multipole 
moments. These moments are analyzed for their far field contributions, which 
are then time averaged to give the power in gravitational waves radiated by 
the source. 

3.1 Retarded Solution to Field Equations 

Let us first start with the equation for the retarded solution to Einstein's 
linearized field equations. 




(3.1) 



If we make the approximation for x » x' then we get 




(3.2) 
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This approximation is valid for a source with any velocity. 

To facilitate the expansion we switch to the TT gauge and find the Fourier 
components of h: 

hj^i^^ d'^' J ^M^^ k)e-(*-^^)e'>^-'. (3.3) 

Next we expand the exponential. 

One can equivalently expand the stress energy tensor around the retarded 
time tr — t — r, 

Tki{^',t-r+'^) « ni{t-r,x') (3.4) 
r 

+x'wd,^T,i + T,, + . . . 

We then collect the terms around the powers of the coordinates into the 
multiple moments. 



3.2 Mass Multipole Moments 

Gravitational radiation in General Relativity is only quadrupolar and higher, 
so dipole or monopoles terms do not create gravitational waves. This is due 
to the photon being spin-2 and massless, thus the only spin states of j = 2 
are allowed. This also reveals itself in the first order solutions of the field 
equations, although it is more general than that. In this first section we write 
down some of the moments from the gravitational field, given by: 



M 



J SxT^ monopole, (3.5) 
= j d^xT^^x^ dipole, (3.6) 
M^'' = J d^xT^^x^x'' quadrupole. (3.7) 

From the equation of conservation of mass/energy in Minkowski space for 
linearized gravity we get, 

a^T"^ = 0, (3.8) 



31 



leads us to the additional restrictions: 

M = (3.9) 

= P^' = J d^xT^^ momentum (3.10) 
= P^ = (3.11) 

M^'^ = 2 / (3.13) 



Of importance arc the higher order derivates of the quadrupolar terms, as 
they lead to gravitational waves. 



3.3 Solutions for h^^y 

If we expand the solution of h^j^y in terms of the retarded time we find that 

2(7 

/loo = — (M-n,P^), (3.14) 

hoj = ^P,-, (3.15) 

h,k = —[-mk{M-n,P') + f,ki^Q'"''). (3.16) 
where the quadrupole tensor is given by, 

Qjk = M,fe - i,7,feM„", (3.17) 

Qjk = y rf^a;Too {^jXk - ^^jkr'^^ , (3-18) 
and the projection tensor is given by, 

fjklm — ^jl^km — ^^jk^lm, (3.19) 

n^-fe = Vjk-njUk, (3.20) 

n, = X (3.21) 
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For our transverse-traceless (TT) gauge we take n = (0,0,1), and our 
equations simplify to, 

hjk = ^ Qjk. (3.22) 

Nice. Here we have ignored the monopole and dipole terms, as they don't 
contribute to the gravitational radiation. 



3.4 Power Radiated 

Now for the power radiated by a source: we must deal with both the energy 
from the source and energy in the waves. First we note that the effective 
stress energy tensor is a sum of both the gravitational wave contribution and 
the source, 

V = (Tf^u + til;), (3.23) 

1 



SttG 



(^R,u - ^gf^uR^ , (3.24) 



where g is the determinant of the spacetimc metric. If the second equation is 
solved using the Christoffel symbols and assuming a fiat background metric 
we find the gravitational wave portion, t^j^, (known as the Landau-Lifshitz 
pseudo stress energy tensor) is given by. 



tii; = 3^ d.h^f^d^h'^^ (3.25) 
Averaging this gives the effective tensor, also known as the Isaacson tensor: 

= (tiu), (3.26) 



c 



,4 



327rG 



d^Kfsduh'^") . (3.27) 



Conservation of the stress-energy tensor is given by the covariant deriva- 
tive, 

D.T^';^ = 0, (3.28) 

which in a flat background is 

doT^^ + d.T^l ^ (3.29) 
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when integrated over a volume gives, 

/ d'x {d,T^l + d.Tll) = 0. (3.30) 
The 00 component is the energy density, so we find, 

L = y d^xd.T^l (3.31) 
when evaluated on the surface gives. 




(3.32) 



This is a statement of Green's theorem. 
Prom this we find the fiux is given by, 

Fg^ = cT™, (3.33) 

where the speed of the wave is c. 

The average gravitational wave power radiated from a source is given by. 

Law = j d^x T^^n\ (3.34) 

The final result is given by, 

L,^ = ^ {Q,k&') , (3.35) 

which was first derived by Einstein a long time ago. In this case Too has 
units of mass per unit volume, rather than energy per unit volume. Note c 
has been reinserted. 



3.4.1 Plane Waves 

Assuming the source very far from the detector, we can make a plane wave 
approximation. For plane waves in the TT gauge the energy flux can be 
calculated to be, 

Fgw = -^fhlms^ (3.36) 

hrms = ^\ {h\ + hl). (3.37) 
34 



From this /i, which is hj-msi can be determined if the luminosity and the 
distance to the source are known, 

h = (3.38) 
Trr/ 

Again, this is vahd if the source of the waves is distant enough to approximate 
a plane wave solution. 



3.5 Astrophysical Sources and the Stochastic 
Background 

3.5.1 Binary Systems 

Let us consider a binary system of equal mass m in circular orbit of radius 
R, with the masses separated by a distance of 2R. The moments of inertia 
can be easily calculated and the time derivatives found using the equation 
above, which leaves us with: 

Law = ^m^i?^a;^ (3.39) 

= 4173^ (^^)^' (3-40) 

for the last equality we have used u^R = Gm/{2R)'^. If we seperate our 
factors of G and c we get, 

32c5 fGmuV^^^ 



= 3.021 X lO'^'W — — , (3.42) 
where T is the period in hours. 



3.5.2 Stochastic Background 

There are a number of sources of a stationary, isotropic gravitational wave 
background, to include: 

• Unresolved whited dwarf and neutron star binaries. 
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• Supernovae. 

• Quantum fluctuations from post-inflation relieating. 

• First order cosmological phase transitions and subsequent turbulence 
from bubble interactions. 

• Topological defects: cosmic strings. 

The universe is expected to be rich in a number of these sources, with the 
possibility of cosmic strings as well. Cosmic string loops are a very discernible 
source due the very flat spectrum in frequency space. The one-scale model 
ensures that the radiation from past epochs is added in such a way that the 
spectrum is very flat. Current loops have the effect of creating a peak at 
approximately the inverse Hubble time of decay. 



The stochastic gravitational wave background can be characterized by the 
dimensionless quantity, 

n,M) - (3.43) 

where = 3H^/{8itG). 

Now we relate to ^Igw The equation for the gravitational wave energy 
density from the waves is. 

From eqn. 12.231 we find, 

= TTTT^ / d{\nf) {27: ff hlif). (3.45) 



If we take the derivative with respect to In/, 

i^nffhiif). (3.46) 



din/ 1671^ 
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Finally we divide by the critical density, 



1 dpg^{f) 

Pc dlnf 

^gui (/) 




1 



■c 



(3.47) 



(3.48) 



fS,{f). 



(3.49) 



3.6 Millisecond Pulsars 



An ingenious method for detecting a stochstic background of gravitational 
waves is the use of millisecond pulsars [53]. The frequency of millisecond 
pulsars is measured for extended periods of time, and all known effects on 
the frequency are removed. The resulting residual frequencies are analyzed 
for evidence of a stochastic background. 

The low frequencies of this method, on the order of inverse 20 years, are 
the only method of detecting the stochastic background predicted in many 
models. This range is below that of LISA, and puts tangible limits on the 
background of cosmic strings. 

3.7 Alternate Theories of Gravity 

All viable metric theories of gravity have gravitational waves |108j . For 
GR the photon is a massless spin-2 particle and thus there are only the 
plus and cross polarizations. Also, the lowest order contribution is from the 
quadrupole moment while other theories of gravity, to include the scalar- 
tensor theory, include a dipole moment [lOSj . Thus LISA could potentially 
be a test of GR for a binary mass in-spiralling system. 

Other effects include a difference between the speeds of light and gravi- 
tational waves, and up to six possible polarizations aside from the aforemen- 
tioned two. 

What follows is a brief summary of dipole radiation taken from Will, 1993 [108j . 
If we define the dipole moment of the self-gravitational binding energy of 
two bodies as. 




(3.50) 



a 
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where is the self gravitational binding energy of object a given by, 

fP^^^d'xd'x'. (3.51) 
So at least representatively one can write, 

Lew = l^^D (d • d) , (3.52) 

where is a constant that depends upon the theory. 

Scalar-tensor theory: As an example we look at a binary system with 
Brans-Dickc (BD) scalar-tensor theory. For the luminosity we have the dipole 
contribution, 

Lf^ = l^D^': (3.53) 

where Q = Qi/rrii — 0,2/ m2 is the difference in the self-gravitational binding 
energy of two bodies. For BD theory ko = 2/(2 -\- u), which is very small 
given the current solar system bound with uj > 500. 

Even though this value is quite small given the current limits, it is still 
discernible as an additional source of gravitational wave energy. 
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Chapter 4 

Gravitational Wave Detectors 



4.1 Gravitational Waves on the Detector 

A point which the reader should take note of is the extremely small amplitude 
of the signal produced by gravitational waves. One way to understand this is 
via the Einstein equation Gfj_y = IQnGT^i, which leads to Eqs. 12.191 and 12.71 
Eq. 12.191 indicates h Ax/x is a strain, and is often referred to as the 
"dimensionless strain of space" or the "time-integrated shear of space" [13]. 
Eq. 12. 71 shows that T^'^ is comparable to the stress. If we use the analog of the 
stress-strain relationship of a solid we find the Young's modulus is given by 
c'^/(167rG); this is a huge number! In other words, spacetime is very "stiff" 
and strongly resists the stress produced by a passing gravitational wave. 

Thus the signal on any detector is going to be minuscule, even for a very 
luninous source. In order to extract this signal the sensitivity of the device 
must, apparently, be very high. This sensitivity then leads to problems with 
background noise, especially on the earth. A typical strain is /i ~ 10~^^, 
which for a detector with arms of length 10 m is a sensitivity on the order of 
a nucleus. For LISA with arms of length 10^ m the length of the detectable 
oscillations is the size of an atom. Getting a signal from the noise at these 
small values is indeed no small feat. 

In general there are two different types of detectors: resonant-mass de- 
tectors and interferometers. Resonant-mass detectors were first pioneered by 
Joseph Weber in the 1960's. Today the detectors are four orders of magni- 
tutde more sensitive than Weber's original bars, but even this sensitivity only 
allows the detection of very powerful emitters within the galaxy and local 
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galactic neighborhood. Given the relative rarity of such events, it is unlikely 
they are to be detected. On the positive side, resonant-mass detectors can 
be fabricated on small scales when compared to the interferometers. 

Interferometers were first theorized for use in 1962 [70] but their complex- 
ity and cost prevented implementation for the next thirty years. The current 
detectors (LIGO, VIRGO, GEO600, and TAMA) are all taking data now, 
but have had no detections. These are all ground based detectors and must 
deal with noise from the earth, especially at frequencies from 1-10 Hz. This 
means that ground based detectors are limited to frequencies above 10 Hz. 
The first space based detector, LISA, is planned to launch early in the next 
decade. Being free from the earth's seismic noise, this detector is planned to 
have a detection band in the millisecond range. 

For details on the calculations of strain from single and multiple detectors 
see Appendix C. 

4.2 Ground Based Detectors 

Ground based detectors typically come in two varieties: resonant bars and 
laser interferometers. Their detection bands are ultimately limited by seismic 
noise from the earth, so they typically cannot go below 10 Hz. 

First we discuss the resonant-mass detectors. These are of less interest 
in our study as their frequency range is of the order 700-900 Hz [22], which 
corresponds to very small loops. This corresponds to loop sizes for which our 
simplified model is not necessarily accurate. 

The laser interfermometers have a much larger sensitivity range, but are 
still limited by seismic events. This puts their range fairly high. 

4.3 Laser Interferometer Space Antenna 

The Laser Interferometer Space Antenna LISA, is in the final stages of pre- 
fiight and is scheduled to launch in the near future. It is a technological 
marvel, using three satellites in an equilateral triangle with sides of 5 x 10^ 
km . It is to be in orbit around the sun at the earth's radius 20° behind the 
earth, or a distnce of 5 x 10^ km . The triangle of satellites are dipped an 
angle of 60° from the ecliptic [221 [63l |Ml EOl [93] • 

The length of the arms correspond to a peak sensitivity of about 10 mHz, 
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see for example 15.11 The frequencies below about 0.1 mHz are determined 
by the drift of the spacecraft, which every attempt is made to minimize. 
Our plots show theoretical values of the sensitivity for those low frequencies, 
although it is by no means certain they will be attained. 

The range of frequences in LISA's detection band is quite rich. A number 
of astrophysical sources to include emission by white dwarf binaries and the 
inspiral of black holes and stars with supermassive black holes can be found 
here. Also included are the most luminous events in the universe: the inspiral 
of two supermassive black holes. 

Some of the signals are so rich that they form a signal that washes out 
all other signals. This "confusion noise" is expected in the higher frequency 
range from white dwarf binaries in the galaxy. Several studies have been done 
on this [76], |37] and their results are indicated on our plots of sensitivity. 
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Chapter 5 



Gravitational Radiation from 
Individual Cosmic Strings 

Perhaps the most exciting aspect of gravitational waves from cosmic strings 
is that the radiation from the loops is a series of harmonic modes. This allows 
even small signals to be extracted from the background. Also, for the light 
strings, the quality factor is very large which ensures the frequencies of the 
harmonic modes stay constant during detection. Given that the maximum 
time frame for detection is on the order of three years, we can safely assume 
the signal is stationary in time. 

There has been much investigation into the radiation from cusps, kinks, 
and at reconnection points [51]. These sources are typically directed bursts 
of gravitational waves and are not harmonic in nature. They also tend to be 
at much higher frequencies than the fundamental modes of the larger string 
loops. For these reasons they are not included in this study. 

In this chapter we first re-analyze the equations of motion for a cosmic 
string loop and determine that the motion is indeed a series of harmonic 
modes. Subsequently the distribution of the loops is estimated using the 
distribution of dark matter as a model. It is assumed that the cosmic string 
loops follow the dark matter as they are bound gravitationally to the galaxy. 
Finally the power and spectra of an individual loop are estimated using 
numerical calculations. 
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5.1 Linearized Gravity 

Let us first start with a review of tlie equation for the retarded solution to 
Einstein's linearized field equations, 

J I-X -X I 

For a cosmic string we get, 

rpt^u ^ f dadT^-i'^^dAX^'dBX^5^^\x - X)), (5.2) 

J 

where X'^(ct, r) = = {T,r{r,a)), and = g^^dAX'^dBX" . Substituting 
values and assuming a locally Minkowski spaceime we get: 

T^. = fij da{N^V, - iV^iV,)5(3)(x' _ r{t',a')), (5.3) 

where V^" = (1, r^), N^' = (0, r^), and t' = r. 

If we make the approximation for x >> x' then we get 

V(x,t) = — / rfVf^,(x',t-r + ). (5.4) 

r J r 

We then use Eq. 15.31 to find our /i^j, at a given direction and distance from 
the source. 

In principle this equation can be used to find the field far from a cosmic 
string loop. In practice, the waves from a cosmic string loop are harmonic 
so we can simply use plane wave solutions and sum over several wavelengths 
to find the luminosity. This is detailed in later sections of this chapter. 

5.2 Periodicity of Cosmic String Loop Mo- 
tion 

We reiterate previous results from Chapter 1 and then add a final calcula- 
tion which involves cosmic string loops. This last part is important as it 
characterizes the loop motion as a series of harmonic modes. 
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Cosmic string motion can be characterized by paramatcrizing the mo- 
tion on a two-dimensional surface called a worldsheet. With the coordinate 
parameterization given by: 

x^" ^x^ia"), a = 0,1, (5.5) 

where the a;^ are the four spacetime coordinates and the a" are the coordi- 
nates on the worldsheet of the string. In general (7° is considered timehke 
and (7^ is considered spacelike. On the worldsheet the spacetime distance is 
given by, 

ds"^ = g^iydz'^'dx" , (5.6) 

= ^abda'^da', (5.8) 

where •jab is the induced metric on the worldsheet. 

The equations of motion for the string loops can be found via the Nambu- 
Goto action: 

S^-fx j y/^d'^a, (5.9) 

where 7 is the determinant of 7a(,. Minimizing the action leads to the follow- 
ing equations, 

where F are the Christoffel symbols. 

In flat spacetime with g^^, = rj^i, and F^^ = the string equations become, 

^ ' ^7"'S;')=0. (5.11) 



Now we use the fact that the Nambu action is invariant under general 
coordinate reparametcrizations, and we choose a specific gauge. This does 
not disturb the generality of the solutions we find, since the action is invariant 
under these gauge transformations. Here we choose: 

701 = 0, (5.12) 
700 + 711 = 0, (5.13) 
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this called the "conformal" gauge, because the worldsheet metric becomes 
conformally flat. 

Our final choice is to set the coordinates by 

i = = (7° (5.14) 
(7 = (7\ (5.15) 

so that we have fixed the time direction to match both our four dimen- 
sional spacetime and the worldsheet. This is sometimes called the "temporal" 
gauge. This leads to the equations of motion, 

i-x' = 0, (5.16) 
x^ + x'^ = 1, (5.17) 
x-x" = 0. (5.18) 

The first equation states the orthogonality of the motion in these coordi- 
nates, which leads to x being the physical velocity of the string. The second 
equation is statment of causahty as the total speed of the motion is c, which 
can also be interpreted as, 

da = -H= = d£;//., (5.19) 
V 1 — a;^ 

E = 1^1 ^^== ^l^jda. (5.20) 

So the coordinates a measure the energy of the string. 

The final equation is a three dimensional wave equation. This equation 
determines the behavior of the strings and has a general solution of the form: 

x(a,t) = i[a(a-t) + b(a-t)], (5.21) 

where constraints from the equations of motion give, 

a' = b' = 1. (5.22) 

We can also write down the stress-energy tensor, momentum, and angular 
momentum in this choice of gauge: 

T'^" = /X J dalx^x" - a;V]5(=^)(x - x(cT,t)), (5.23) 
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which leads tot he energy given above, as well as, 

p = // y da±{a,t), (5.24) 
3 = jj, da x{a,t) X x{a,t). (5.25) 



It is found for loops that the large angular momentum tends to reduce rota- 
tion rates dramatically. 

5.2.1 Loop Motion 

For a loop of length L we have the identity, x((T + L,t) — x((j, t). Thus the 
functions from our general solution also have the same periodicity: 

a{a±t + L) = a{a±t), (5.26) 

and similarly for h{a, t). Prom above, the smallest value of t that is periodic 
is t — L/2 since, 

a([(T + L/2] ±[t + L/2]) = a((7 ± t). (5.27) 
Hence the fundamental frequency is given by 

/^i^f. (5.28) 

where c has been reinserted for clarity. 

The general solution may be represented in terms of its Fourier compo- 
nents as, 

/ N P 

x(cr, t) = Xo H 1 

A* 

oo 

+ ^ (A^e^^'^'^l'"*]/-^ + B^e-*2'^"[''±1/-^) , (5.29) 

n=l 

where a± — a ± t. So in general, one can write the string motion in terms 
of notmal modes which arc harmonic in time. This fact is important and 
dramatically increases the detectability of the normal mode oscillations of 
the string loops. 
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5.3 Distrubtion of Loops 



The average number density of loops of size L is n{L), computed numerically. 
Thus the total number density of loops is given by, 

= j n{L)dL (5.30) 

and the average distance between loops is {N /V)^^^^ . 

The loops are assumed to be clustered around galaxies within the dark 
matter halo. From the average distance between galaxies the total number 
density, the number of loops within the galactic halo is estimated. 

The distribution detected from the Earth is a function of both the length 
of loops and their distance from the earth. The number density in the galaxy 
will be matched to the dark matter halo, given by the NFW density distri- 
bution pNFwir) IZ5]. The number density is then a function of the length of 
loops and distance from the center of the galaxy n{L,r). 



5.3.1 Calculations of Loop Density 

We start by figuring out the size of the loops at the fundamental frequency: 

L = 2c/ h (5.31) 

We then find the time in the past when loops formed and have decayed via 
graviatational waves to the size now given by 2c/ fi. Using the one-scale 
model loops form at a given time with size L acH~^{t) and start decaying 
at a rate L = —•jGfi. So loops created at a time tc will be of the size L(tc, t) 
at the time t given by the equation, 

L{tc, t) = acH-\Q - -fGfx{t - Q/c, (5.32) 

From the present time to we insert the value of L to get: 

2c/ /i = acH-\tc) --fGfx{t-tc)/c. (5.33) 

The number density as a function of the time of observation, to, and time 
of creation, tc within a At in log scale of 0.1, is given by: 

„it.x>-^(!^rm)\ (5.34) 

a \ c J \a[to) J 
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Using the time tc we solve numerically to find the number density and taking 
we find the number within the galaxy. 
If we assume the average distance between galaxies is 5 Mpc, the number 
of loops within the galaxy for Gfi = 10~^^ and a fundamental frequency of 
1 mHz is approximately 10^. This number is distributed within the halo of 
each galaxy, which leads to an overall increase within the halo and a large 
increase in the density near the center of a galaxy. 



5.3.2 Number Density in the Milky Way 

The density of galactic string loops is matched to the dark matter distribution 
in the Milky Way using the NFW density PNFwif^)'- 

pNFw{r) = (5.35) 
a;(l + xy 

X = (5.36) 

where po and rg are determined by observation. Representative values are 
given by [59] using the favored model with rs=21.5 kpc. From above we see 
that for r << r^, p oc and for r >> r^, p oc r~^. 

The NFW mass density can be matched to the number density of the 
same form: 

where again x = r/rg. The number density constant is matched to this 
curve by integrating over the halo volume up to a truncation value of rt and 
normalizing to the total number of string loops A^: 

rn 

N{L) = / n{r,L)d^x, (5.38) 
Jo 

, ^ A^(^) (n + , , 

ng{V) = p ^ ^ ^ * 7 \ -■ (5-39) 

47rr? 



s 



in + r.)ln (^) - n 



In general is of the order 100 kpc. The value of ns{L) must be determined 
for each length (or equivalently fundamental frequency) of loop, in log bins 
of 0.1. 
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Sample Calculation of Number Density Distribution 

A reprentative value from the previous example of strings of tension G/i — 
10~^^, N — 10^, rt/vs — 10, and = 21.5 kpc, at the earth's distance form 

the galactic center, r = 8 kpc, we find n ~ 9.4 x 10~^ kpc~^. Inverting and 
taking the cube root we find n~^^^ ~ 2.2 kpc. Thus on average the distace 
to the nearest loop is appoximately 2.2 kpc. 



5.3.3 Distribution of Loops Near the Earth 

The signal at the earth will be a function of both distance from the earth 
r' and size L. The number of loops dN of length L a distance r' from the 
Earth and R from the center of the milky way, within an interval dr' and a 
solid angle dQ' is given by: 



where 



dN{L) = n{L, R + r')d'r\ (5.40) 



n{L, R + r') = n{L, r', 9', cj>') = ,!\'^_^\,, , (5.41) 



and 



, _ + i?2 _ 2Rr' sin 9' cos (f)' 
rs 

R is the distance of the Earth from the galactic center and is defined above. 
Here 9' and 0' define the angle measured from the earth. Here 0' = is the 
direction of the galactic center and 9' = 7i/2 defines the plane of the galactic 
disk. The signal in the direction of the galactic center has an increased 
luminosity over other directions, especially directly away from the center. 
The number of particles per unit solid angle from the earth is defined by, 

dN n 

This is useful for determining the direction of highest gravitational wave 
luminosity. 
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5.3.4 Results for Number Density and Distributions 

The number density of loops in the Milky Way is strongly influenced by the 
size and mass density of the cosmic string loops. Tables I5.1H5.3I give results 
for the number of loops in the Milky Way for various parameterizations of 
the string loops. The results clearly show the large predominance of loops 
for very light and large strings. 



5.4 Power from each loop 

A loop will radiate power in each mode n modeled by, 

En{r') (X ^, (5.44) 

where £"„ is power radiated per mode from a loop. The total power is 

E = ^Gix' = Y,PnG^l^ (5.45) 

where 7 is found to be between 50 and 100, and the P„ are power coefficents 
of the individual modes. Note that the power radiated is independent of the 
loop size. To first order we ignore directionality of the emitted gravitational 
radiation. 

Each loop has a spectrum given by each of the power coefficients, P„, 
where the sum ^ P„ = 7 is fixed by the fundamental properties of the strings. 
This is an average over a wide array of loops: it is expected individual loops 
vary somewhat from this. Regardless is expected that the most power is in 
the fundamental mode and the lowest modes, while the higher modes have 
significantly reduced power. 



5.4.1 Quality Factor 

For a given mode the quality factor is given by. 



Q. = ^, (5.46) 
= ^- (5.47) 
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Upon inserting values for a given string loop, we find 



Q 



n 



27TUE 
En 

2n fiL 



(5.48) 



(5.49) 



2' 



which simplifies to the result: 



Q 



Ann 



(5.50) 



n 



For light strings with G/i = 10 



we find, 



Q 



n 




(5.51) 



For the fundamental mode P„ is of order unity at most, so Q is large, and 
increases with higher n. Thus the full width at half maximum should be small 
at each mode, and justifies the use of delta functions in other applications. 

Given the extremely large value of Qn it can be assumed that the signal 
from each mode in the loop is monochromatic for the duration of observation 
with regard to decay, with is at most three years. Potential frequency shifts 
due to gravitational lensing, gravitaional doppler effect, and relative motion 
have not been analyzed [31]. 



5.5 Strain Produced by Gravitational Radia- 
tion 



The relative weakness of gravitational waves from cosmic strings allows the 
use of linearized gravity. The stress energy tensor can then be given by the 
Isaacson tensor: 




(5.52) 



and the total luninosity in gravitational waves is given by 




(5.53) 
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Assuming plane wave solutions, we find for the energy flux: 



4.Q 



(5.54) 



Inserting the radiated power Lg^ = En into the equation above and solv- 
ing for hn we find, 



Inserting c we find : 



■c G/i 1 



^ n 9 r 1 

K = 3.095 X lQ-^'^^/p'n [G//(c = 1)] X 



(5.55) 
(5.56) 

(5.57) 
(5.58) 



For Gil = 10 /=! mHz, and r—1 kpc, we find hi ^ 10 within the 
detection limits of LISA. 



5.6 Spectrum of a Single Loop 

Using previous results for the strain produced by a gravitational wave, we 
can add the harmonic time dependence to find the strain at each mode as 
follows: 



where. 



Thus we have. 



Jn V TTC^ 



^n(t) = ^^-e-^^-^"* (5.61) 
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Again, the P„ are the "power coefficients" and has been estimated to go 
asymptotically as n~^^^. From previous we recall that the sum of our co- 
efficients is given by ^ P„ = 7- In practice most of the power is in the 
fundamental mode, so the the power from lower modes drops off, allowing 
truncation of the sum at a reasonable value of n. As an example, we set 

Pn = Pin-^'\ (5.62) 

and Pi = 18. At n = 20 we find ^ P„ ^ 45, giving us a reasonable estimate. 
For a single loop the total strain will be a sum of all the modes, 

oo 

h{t) = J2hn{t), (5.63) 

n=l 

This sum is typically taken to n = 10 to facilitate computational efficiency. 
It is assumed most of the power is in or near the fundamental mode, with 
power dropping off as as mentioned previously. We expect to lose little 
by ignoring the higher modes given their small contribution. Other research 
has been done on bursts and cusps which are effectively part of this high 
frequency tail. 



5.6.1 Estimate of Strain compared to Stochastic Back- 
ground 

For comparison we investigate the strain produced by the stochastic back- 
ground. Of use is the value flgwif) which is given by: 

^""j'hl (5.65) 



where pc is the critical density of the universe. Ho is the current Hubble 
constant, and hf is essentially the Fourier transform of the gravitational 
wave signal . 

From above we get. 



■lHz3/2 



2.239 X lO-^V^s-^'-pT^. (^•^'^^ 
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where h = Ho/100 ^ 74. 

From previous numerical calculations we find Qh"^ = 10~^'^ for f=l mHz, 
a = 0.1, Gfi = 10"^^, and 7 = 50, gives hf ^ 10~^^. Thus the background 
appears to overwhelm the loops farther than 1 kpc from the earth. For 
distances from the earth to a single loop less than about 100 pc, the gravi- 
tational wave signal from the individual loops is greater than the stochastic 
background. The periodic signal from the individual loops makes it likely 
that even strains less than the background can be picked up by filtering out 
the periodic signal from the background. 



5.6.2 Fourier Transform 

The frequency spectrum from a gravitational wave on a detector is found by 
taking the Fourier transform of the detected time signal. In practice this is 
done over a finite time, thus introducing resolution issues even though, in 
principle, the spectrum of a harmonic signal should be delta functions. For 
the n^^ component of the strain we find. 



hnif) 



1 



T/2 



T/2 



^2 ^ 



VTnfn 



T/2 



^^g.2.(/-/„)t 



T/2 



(5.68) 
(5.69) 



which can be readily integrated to yield, up to a phase. 



hnif) 



C\fRn Gjl 1 



sin[7r(/ - /„)T] 



<f-fn) 

fGjx\ 1 sin[7r(/ - fn)T] 
^2 



r fnif - fn) 



(5.70) 
(5.71) 



Note the introduction of the factor of the square root of the integration time 
T. This is done to match our results to the sensitivity curves of LISA [631IM] . 
This leads to the the result, 



Sh{f) = Hf) 



(5.72) 



where Sh{f) is the spectral density, whose square root is plotted on the 
sensitivity curves for LISA. 
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Thus the spectrum from a loop is the sum of the Fourier components for 
each mode, given by, 

oo 



n=l 



c fGfx\ 1 ^ V^sin[7r(/ - /„)r] 



V^TT^ V ^ fnif - fn) 



( G[i\ 1 ^ sin[7r(/ - n/i)T] 



Given that gravitational wave detectors sample the data as a discrete 
time series, we use a discrete-time Fourier transform, as discussed in the 
next section. 



5.6.3 Signal to Noise 

A gravitational wave of Fourier transform incident on a dector has a 
signal to noise ratio of, 

(5iVi?)2 = 4jr°°d/M:^. (5.76) 
This is an upper limit assuming optimal filtering. 

5.6.4 Calculation of Single Loop Spectrum using Dis- 
crete Fourier Transform 

Since the signal detected by LISA will be as a discrete set of points, it is 
more appropriate to take a discrete Fourier transform of a sample signal. 
In general LISA has a sampling rate of /samp = lOHz, which leads to large 
numbers of data points for integration times of one to three years. This large 
number of data points also leads to a very clear signal due to the harmonic 
nature of the gravitaional waves produced by the cosmic strings. 

Given the signal h{j) where j is an integer and the data point cor- 
responding to hit) — h{j/ST), we define the discrete-time Fourier transform 
as. 
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is the total number of data points and k is the integer value in Fourier 
frequency space. Note again the introduction of the square root of the inte- 
gration time T ClS db normalization. 

The gravitational wave signal from the cosmic strings has the form, 

hU) = E sin(27rj^„/iV), (5.78) 

where a;„ is not the physical frequency, but the discrete time frequency, which 
is different than the physical frequency given by /„ = 2nc/L, and P„ is the 
magnitude of the n*^ power mode. Using Pn = Pin~^^^ we find, 

•' n=l 

which is then inserted into 15.771 



5.7 Results 

In this study an estimate of the distances to the loops nearest the solar 
system are calculatd for each frequency and string tension, and from this we 
calculate the signal of the loop from our location, the solar system. 

5.7.1 General Results 

Cosmic string loops have a distinctive spectrum of harmonic modes labelled 
by n, with a high frequency tail that goes as n~^^^. 

We use the one year results for the LISA sensitivity plots and most of our 
simulations are run for one year, aside from the higher frequency, /i > 10~^ 
Hz, loops. In these cases we use one month sampling time. 

Investigations of the parameter space of the loops indicates the heaviest 
detectable loops individually are of tension Gfi > 10~^° for a = 0.1, due to 
their reduced numbers in the Milky Way, see Table 17. 2[ 

For smaller loops the detectable range decreases due to a smaller number 
in the galaxy. Shown in Fig. 15. H is the largest likely signal from an individual 
loop at the frequency 10~^. This signal stands above the noise from binaries. 

For the large loops, a = 0.1, it is found that the smallest loops likely 
detectable are of tension Gfi = 10~^^. Fig. 15.31 gives an example within the 
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log(SA(f)"2 Hz-"-) 




Figure 5.1: Plot of the cosmic string loop spectrum for G/i = 10"^^, a=0.1, 
at a distance of r = 7 kpc. Sampling time is T = 1 year, the sampling rate 
is 0.1 Hz, Pi = 18, and the fundamental frequemcy is /i = 10~^ Hz. The 
first 10 modes are shown. This signal stands above the confusion noise from 
galactic binaries which is shown on the graph. The LISA sensitivity curve is 
given by the thick gray line. 
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binary confusion noise, but the loops that have decayed are detectable at 
higher frequencies. 

5.7.2 Varying Loop Size 

The size of the loops formed as a fraction of the horizon, a, affects the number 
of loops currently radiating. For large a the number of loops remaining is 
large, thus at a given frequency the distance to the nearest loop will be 
shorter. This is indicated in Fig. 15.21 Note the amplitude varies as r~^, 
which shows only a small difference on our log plots. 

The total number of loops is significantly larger for large a. This indicates 
that the total signal from all loops in the galaxy at a particular frequency 
is significantly enhanced for large loops campared to that for the individual 
loops. 

Shown in Fig. 15.21 is a plot of the two frequency spectra with varying a: 
one is 0.1 the other is 10~^. The sampling time is T = 1 year, sampling 
rate is 0.1 Hz, G/i = 10~^^, Pi = 18, r = 2.2kpc, and the fundamental 
frequency /i = 1 mHz. In this case the string tensions differe by four orders 
of magnitude, while the distances differ by a factor of approximately 20. 
The ratio of approximately three orders of magnitude matches the overall 
difference between the signals from the loops. 

5.7.3 Varying String Tension 

The increase in the radiated power of the heavier string loops results in great- 
est variation in signals potentially detectable. Fig. 15.31 shows the difference 
in signal via a string of tension G/i = 10~^^. The heavier string loops are, on 
average, farther away but their larger output makes up for increased distance. 

This is, on average, a typical result. The upper limit occurs at approxi- 
mately G/i > 10~^°, where the number of loops per log frequency bin becomes 
of order one (and smaller). At these large string tension we expect the prob- 
ability of there being a loop within the galaxy to vary from somewhat likely 
to very unlikely. 



58 



-14 r 



-16 - 









I 








-1^ 


















[ -20 





-6 -5 -4 -3 -2 -1 

Figure 5.2: Plot of two cosmic string loop spectra with Gfi = 10~^^: the 
"thin" spectrum a=0.1 at a distance r = 2.2 kpc from the solar system, and 
the "thick" spectrum a = 10~^ with r=8.9 kpc. Sampling time is T = 1 
year, the sampling rate is 0.1 Hz, and Pi = 18. The fundamental frequemcy 
is /i = 1 mHz and the first 10 modes are shown as is the confusion noise 
from galactic binaries. 
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Figure 5.3: Plot of the Fourier spectrum of cosmic string loops with tension 
G/i = 10~^^ and a distance of r = 0.066 kpc. Sampling time is T = 1 year, 
and the sampling rate is 0.1 Hz, a = 0.1, and Pi = 18. The fundamental 
frequemcy is /i = 1 mHz and the first 10 modes are shown. Note the 
heavier loops (G/i = 10~^^) have a much larger signal, in spite of their greater 
distance from the solar system. The confusion noise from galactic binaries is 
also shown. At this fundamental frequency, the loop is not detectable, but 
shifter to higher frequencies it is. 
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5.8 Discussion 



These results dramatically increase the detectable regime of cosmic string 
loops. Due to the constraints by recent observations of the cosmic microwave 
background and millisecond pulsars, light strings appear to be the most likely 
candidates remaining for cosmic string phenomena. These light strings have 
a very low stochastic background signal, which is reduced further if the loops 
fracture into tiny sizes. 

These results indicate that very light cosmic strings are detectable due to 
the number density of loops and subsequent clustering around dark matter, as 
well as the harmonic nature of the cosmic string gravitational wave emissions. 
The signal of an individual loop is a series of harmonic modes with a very 
large quality factor. Thus the frequency spectrum is static over the several 
years of observation, further increasing detectability of the signal. 

This study has taken a somewhat simplified approach to the distribution 
of the loops, simply matching them to the dark matter halo of the Milky 
Way. No attempt was made to take into account redshift or any large scale 
anisotropics in cosmic string production. Neither have we taken into account 
directionality in the gravitational radiation from the string. There has been 
much work in the literature on this subject for particular string configura- 
tions [M]. 

The following chapter gives the results of numerical calculations of the 
spectrum of all loops within the Milky Way. 
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Table 5.1: Number of Loops in the Milky Way for the given fundamental 
frequencies and a— 0.1 



GjJL 


/i (Hz) 10-4 


10-3 


10-2 


10-1 


10-11 


338 


338 


338 


338 


10-12 


1.05 xlO^ 


1.05 xlO'i 


1.05 xlO'i 


1.05 xlO-i 


10-13 


3.29 xl05 


3.33 xl05 


3.34 xlO^ 


3.34 xl05 


10-14 


9.22 xlO^ 


1.04 xlO^ 


1.05 xlO^ 


1.05 xlO^ 


10-15 


1.25 xlO^ 


2.92 xlO« 


3.28 xlO« 


3.32 xlO^ 


10-16 


3.21 xlO^ 


3.94 xlO^ 


9.22 xlO^ 


1.04 xlOi° 



Table 5.2: Number of Loops in the Milky Way for the given fundamental 
frequencies, varying a and G/i. = IO-12 



a 


A (Hz) 10-4 


10-3 


10-2 


10-1 


10-1 


1.05 xl04 


1.05 xl04 


1.05 xl04 


1.05 xl04 


10-2 


3400 


3400 


3400 


3400 


10-3 


1100 


1100 


1100 


1100 


10-4 


385 


386 


386 


386 


10-5 


160 


161 


161 


161 


IQ-*' 


93 


93 


93 


93 



Table 5.3: Number of Loops in the Milky Way for the given fundamental 

frequencies, varying a and G/i = lO-i^ 



a 


/i(Hz) 10-4 


10-3 


10-2 


10-1 


10" 


-1 


3.21 xlO^ 


3.94 xlO^ 


9.22 xlO^ 


1.04 xlOi° 


10- 


-2 


1.02 xlO^ 


1.33 xlO^ 


3.30 xlO^ 


3.76 xlO^ 


10- 


-3 


3.21 xlO^ 


3.97 xlO« 


9.32 xlO« 


1.05 xlO^ 


10- 


-4 


1.02 xlO^ 


1.25 xlO^ 


2.92 xlO^ 


3.29 xlO^ 


10- 


-5 


3.26 xlO^ 


3.97 xlO^ 


9.27 xlO^ 


1.04 xlO^ 


10- 


-6 


1.07 xlO^ 


1.27 xlO^ 


2.96 xlO^ 


3.33 xlO^ 
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Chapter 6 

Multiple Loops in the Milky 
Way 

6.1 Multiple Loops 

The previous chapter went into great detail about the spectrum from one 
loop, so this chapter is dedicated to the actual signal the sum of many loops 
in the galaxy produce. Previous studies of binary systems are relevant as 
they provide a model for detection [761 EZ] • The total of all the loop spectra 
produces a background noise that obscures other signals, while a few of the 
closer loops stand out above this noise. The signal is strongly enhanced by 
the harmonic nature of the gravitational waves from the loops. 

Using the number of loops in the galaxy we found earlier, we integrate 
over the whole sky to sum all of the loops. Of significance is the relatively 
large numbers loops present if the progenitor strings are light. This could 
produce a fairly significant signal, especially luminous in the direction of the 
galactic center. 

The plots of all of the spectra in this work show an average signal akin 
to the confusion noise from the white dwarf binaries, as they obscure a fairly 
significant range of frequencies. It is expected that the noise from cosmic 
strings is similar, with variations in its magnitude and frequency based on 
the string parameters. The larger the string tension the greater the potential 
signal, limited by the the diminished numbers in the galaxy. At Gfi > 10~^^ 
it is found that the number of loops in the galaxy drops significantly (see 
Table El]). 
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6.2 Power and Strain from Loops 



We reiterate here that the flux of gravitational wave energy from a loop in 
each mode n is modeled by, 

Frir') oc ^, (6.1) 

where £!„ is power radiated per mode from a loop. The total power is 

£; = 7G'/x2 = J]P„G'/x^ (6.2) 

where 7 is found to be between 50 and 100, and the P„ are power coefficents 
of the individual modes. Each loop has a spectrum given by each of the 
power coefficients, P„, where the sum ^ P„ = 7 is fixed by the fundamental 
properties of the strings. This is an average over a wide array of loops: 
it is expected individual loops vary somewhat from this. Regardless, it is 
expected that the most power is in the fundamental mode and the lowest 
modes, while the higher modes have significantly less power. 

The relative weakness of gravitational waves from cosmic strings allows 
the use of linearized gravity. The stress energy tensor can then be given by 
the Isaacson tensor: ^ 

~ 327rG\ '^''^ r ^^■^> 
and the total luninosity and the flux in gravitational waves are given by 

- / Tl'^nH^x, (6.4) 
Pg^^cToo, (6.5) 

which are found from the conservation of the stress-energy tensor D^^^ — 
in a flat (Minkowski) background. 

Assuming plane wave solutions in the TT gauge, we find for the energy 
flux: 

F,^ = j^fX, (6.6) 

where h is the root-mcan-square of the strain. 
Solving for the rms strain we get, 



>, 2 GFr 
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where each loop has a flux given by, 




(6.8) 



(6.9) 



Thus from the total flux we can find the rms strain at the detector. 



6.3 Total Galactic Flux 



We then sum over all of the loops in the galaxy, grouping the results in 
bins of A/ k, 1/T where T is the total integration time. Typically this is 
between 1 and 3 years, for these calculations we use 1 year. This value is 
indicated on the LISA sensitivity curves, which shows the minimum hrms 
that is detectable for the given integration time. 

6.3.1 GW Signal from Galactic Loops 

Another key ingredient is the gravitational wave signal from loops within 
the galactic halo. These loops are not resolvable on an individual basis, but 
do contribute to the background. The isotropic and stationary background 
from both current and evaporated loops has been calculated previously. To 
this background we add the signal from loops within the dark matter halo of 
the milky way, which is not isotropic from the solar system. The results are 
given in a later chapter. 

So for the m}'^ mode of the flux at a frequency we sum all the contri- 
butions in the galaxy to get the spectrum: 



where n(x, L) is the number density of loops and the primes denote coordi- 
nates centered on the solar system. Both L and are related by L = '^c/ fi. 
We repeat here for clarity. 




(6.10) 



n{Ly,9'A') 



ns{L) 




X 



y/r'^ + i?2 - 2Rr' sin 6' cos 0' 



(6.12) 



s 
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The number density constant is matched to this curve by integrating 
over the halo volume up to a truncation value of (about 150 kpc) and 
normalizing to the total number of string loops A^: 

N{L) = J\{r,L)£x, (6.13) 

ns{L) = p y J. (6.14) 



47rr? 



in + r,)ln - n 



In general rt is of the order 100 kpc. The value of ns{L) must be determined 
for each length (or equivalently fundamental frequency) of loop. 

One can also calculate the expected gravitational wave energy flux from 
a region of the sky: 

= Frn{L,ry'. (6.15) 
Results of numerical calculations of the flux are given in Figs. 16.11 and 16.21 

6.4 Strain Spectra for Light Loops 

From the total flux at a given frequency we find the strain, 



2_ GFgJJn^'''' 

/nV 7rc3 



hifn) = ^\ r^rr ^ (6-16) 



which should be interpreted as the rms strain at a frequency /„. 
Results of numerical calculations are plotted in Figs. 16.31 and 16. 4[ 



6.5 Heavier Loops 

The average spectra calculated above do not show individual string signals 
that may stand out. For the light loops the number of loops is quite large 
so the average is an accurate description of the likely signal. For heavier 
loops Table 15.11 shows the sharp decrease in number of loops at a particular 
frequency. This reduction requires an analysis of each loop and its potential 
signal, while using an average is less accurate. 
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Figure 6.1: Plot of the cosmic string loop gravitational wave flux at the Earth 
for large loops a = 0.1. The uppermost curve is string tension Gfi = 10~^^ 
down to 10"^*^ in increments of 10^. 
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Figure 6.2: Plot of the cosmic string loop gravitational wave flux at the Earth 
for small loops a = 10~^. The uppermost curve is string tension Gfi = 10~^^ 
down to 10"^*^ in increments of 10^. 
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Figure 6.3: Plot of the cosmic string loop strain spectrum for large loops 
a = 0.1 in the galaxy. The top curve is of string tension Gfi = 10~^^ and 
the bottom curve is Gfi = 10~^°, in increments of 10^. Also included are the 
LISA sensitivity curve with an integration time of 1 year, and the galactic 
white dwarf noise. 
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Figure 6.4: Plot of the cosmic string loop strain spectrum for small loops 
a = 10~^ in the galaxy. The top curve is of string tension Gfi = 10~^^ and 
the bottom curve is G/i = 10~^°, in increments of 10^. Also included are the 
LISA sensitivity curve with an integration time of 1 year, and the galactic 
white dwarf noise. 
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Figure 6.5: Plot of the cosmic string loop strain spectrum for small and 
large loops in the galaxy. It is assumed 90% of the loops shatter into small 
loops of a = 10~^. The top curve is of string tension Gfi = 10~^^ and the 
bottom curve is Gfi = 10~^°, in increments of 10^. Also included are the 
LISA sensitivity curve with an integration time of 1 year, and the galactic 
white dwarf noise. 
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For loops of tension Gjj, > 10^^° the number of loops in the galaxy is 
likely of order 10^ overall, and though their distribution is roughly that of 
the dark matter the small numbers require more care in analysis. 

In each frequency bin there are few loops, so its distribution is random 
with probability given by the NFW distribution. Thus the signal may be 
very low-perhaps not detectable-if the loop is far away, or very large if the 
loops is close to the earth. 

6.6 Discussion 

The detectable strain of light loops is strongly enhanced by the periodic 
nature of the signal and the clustering of string loops within the galactic 
halo. This effect is greater for large loops, but even small loops show a 
detectable signal for LISA. 

For large loops, detectable string tensions low as G/i = 10 

three orders of magnitude lower than previous estimates [21]. Even if the 
loops shatter into very small loops, the detectable signal is still appreciably 
larger with string tensions down to Gfi = 10^^^. 
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Chapter 7 



Gravitational Wave 
Background from Cosmic 
Strings 



7.1 Introduction 

In this chapter we calculate the spectrum of background gravitational ra- 
diation from cosmic string loops. The results of these calculations confirm 
the estimates made in an earlier exploratory analysis, that also gives simple 
scahng laws, additional context and background [16] . Here we add a realistic 
concordance cosmology and detailed numerical integration of various effects. 
We estimate that the approximations used to derive the current results are 
reliable to better than about fifty percent, and are therefore useful for com- 
parison of fundamental theory with real data over a wide range of string 
parameters. 

Cosmic strings radiate the bulk of their gravitational radiation from 
loops [H I106[ [99l [Ml 1103] . For a large network of loops we can ignore di- 
rectionality for the bulk of their radiation which then forms a stochastic 
background [1]. Each loop contributes power to the background: 

^ = 7G/.^c, (7.1) 

where 7 is generally given to be on the order of 50 to 100 |l6l 11061 [T71 HOl [2]. 

Gravitational radiation from strings is characterized for our purposes by 
the dimensionless luminosity of gravitational waves 7, which depends on the 
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excitation spectrum of typical radiating modes on string loops. Loop for- 
mation is characterized by a, the typical size of newly formed loops as a 
fraction of the Hubble scale H{t)~^. Estimates based on numerical simula- 
tions have in the past suggested very small values, leading to the hypothesis 
that a ~ Gfi or smaller |89l [12], |6] ; those ideas suggest that all but a frac- 
tion ~ Gfi or smaller of loops shatter into tiny pieces. Recent simulations 
designed to avoid numerical artifacts suggest a radically different picture, 
that in fact a significant fraction of loops land on stable non-self-intersecting 
orbits at a fraction ~ 0.1 of the horizon scale [lOH 11021 [8il [71] . Our study 
is oriented towards this view, which leads to a larger density of loops and a 
more intense background for a given string mass. 

For units of the gravitational energy in a stochastic background a con- 
venient measure we adopt is the conventional dimensionless quantity given 
by 

= (7.2) 
pcdlnf 

where pgw is the energy density in the gravitational waves and pc is the 
critical density. Qgu]{f) gives the energy density of gravitational waves per 
log of the frequency in units of the critical density, flgw is proportional to 
the mean square of the strain on spacetime from a stochastic background of 
gravitational waves. 

There has been much interest in cusps and kinks as sources of gravita- 
tional radiation from cosmic strings [271 [28l [29l [88] . this paper an estimate 
of the radiation power from cusps and kinks is included as higher frequency 
behavior in the loop gravitational wave spectrum, but the beaming and bursts 
of gravitational waves is not discussed explicitly. As estimated in [16], in the 
regime of very light strings, observable bursts are expected to be rare, and 
harder to detect than the stochastic background. 

Finally we would like to state some additional uncertainties associated 
with the model of strings described below: 

• There is the possibility of multiple scales in loop formation: some string 
theories suggest multiple stable intersection points and loops connected 
by these intersections. This may strongly affect the accuracy of the 
one scale model, and indeed some of these theories also have signifi- 
cant other modes of decay that invalidate the entire framework used 
here 1501. 
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• The reconnection probability p can be less than one. This may affect 
loop formation and the gravitational wave background. In general it 
is thought that this effect always increases the predicted background 
although it is not agreed by how much [35|. 

• There may be transient behavior of cosmic strings due to the formation 
process of cosmic quenching. The model below assumes that transient 
effects have long since settled out by the time our model begins; this 
seems likely at LISA frequencies which correspond to initial loop sizes 
much larger than the horizon at the end of inflation 

Further numerical string network studies will be needed to resolve these is- 
sues. 

7.2 Model of String Loop Populations 
7.2.1 Loop Formation 

The behavior of strings on cosmological scales has been thoroughly discussed 
in the literature [M HOl HI 113 [Ml [M |95l El |12]. By the Kibble mech- 
anism strings are created as a random walk with small coherence length, 
and they evolve following the Nambu-Goto action. The strings form in 
large threadlike networks jlOlj . and on scales longer than the Hubble length 
cH{t)~^ the network is frozen and stretches with expansion. The strings 
move with speed c, and they interact to form loops with a probability p 
by different mechanisms: one involves two strings, the other a single string 
forming a loop. These loops break off from the infinite string network and 
become the dominant source of gravitational waves. They oscillate with a 
fundamental frequency given by / = 2c/ L, where L is the length of the loop. 
In general the loops form a discrete set of frequencies with higher modes 
given by fn = 2nc/L. 

The exchange probability p measures the likelihood that two crossing 
string segments interact and form new connections to one another. If p=0 
then the string segments simply pass through each other; but if p=l then they 
can exchange "partners", possibly forming a loop. The value of p depends 
upon the model [351 1121 [36] ; it is close to unity in models most commonly 
discussed but in principle is an independent parameter (and in a broad class of 
models, the only one) that can modulate the amplitude of the spectrum for a 
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given string tension. In this work it is taken to be p=l and is parameterized 
as part of the number of loops formed at a given time [ISl I1U31 1^ . The 
number of loops is normalized to the previous results of R.R. Caldwell and 
B. Allen [18], which is described in detail later in this section. 



7.2.2 Loop Radiation 

For the dynamics of a string the Nambu-Goto action is used: 



S = -fi ^/-g(^)dTda, (7.3) 



where r and a are parameterized coordinates on the world sheet of the loop, 
and the determinant of the induced metric on the worldsheet. The 

energy-momentum tensor of the string can be given by |106j . 

T^'^(t,x)=/i [ dafi'^x" -'^^]5^^\x-x{t,a)), (7.4) 
J \ da da J 

where r is taken along the time direction. The radiated power from a loop 
at a particular frequency, /„, into a solid angle fl is, 

"^^^ 4vrG/„^ fr* (/„, k)T'^'^{fn, k) - 1 |T;(/„, k)A . (7.5) 



Here T^uifn, k) is the Fourier transform of T^uit, x) [106] . Dimensionally the 
power is given by G^^ [Mj with numerical coefficients determined by the 
frequency mode and the luminosity of the loop. 

In general for the power output of a loop at a given mode we can write 
En = PnGjJ'^c |106[ [T8]. where directionality is ignored by effectively inserting 
the loop into a large network of loops and the P„ are dimensionless power 
coefficients. Note that c has been reinserted. The power coefficients are 
generally given the form P„ oc n~^^'^ |106] which will create an f~^^^ high 
frequency tail to the power spectrum of a loop. The total power from a loop 
is a sum of all the modes: 

oo 

E = ^P^G^^c, (7.6a) 



n=l 



-fCfi'c. (7.6b) 
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Thus 7 measures the output of energy and is labeled the "dimensionless 
gravitational wave luminosity" of a loop. For most calculations we use the 
value 7 = 50 [2]. 

For the majority of these calculations we assume only the fundamental 
mode contributes, effectively a delta function for P{f), and / = 2c/ L is the 
only emission frequency of each loop. This assumption is relatively good. To 
check this, some calculations include higher frequency contributions and show 
that the effect on the background is small overall, with almost no difference 
at high frequencies. This is discussed in more detail in Results. The high 
frequency additions also represent the contribution from kinks and cusps, 
which add to the high frequency region with an /^^^^ tail from each loop. 
This high frequency dependence is the same as that for the high frequency 
modes, so their contributions should be similar, particularly for light strings. 
Light strings contribute at higher frequencies because the loops decay more 
slowly and more small loops are able to contribute at their fundamental 
frequency longer throughout cosmic history. 

The populations of loops contributing to the background fall naturally 
into two categories [IB]: the high redshift "H" population of long-decayed 
loops, and the present-day "P" population, where the redshift is of the order 
unity or less. The spectrum from the H population, which have long decayed 
and whose radiation is highly redshifted, is nearly independent of frequency 
for an early universe with a p = p/3 equation of state. The P loops have a 
spectrum of sizes up to the Hubble length; they create a peak in the spectrum 
corresponding approximately to the fundamental mode of the loops whose 
decay time is the Hubble time today, and a high frequency tail mainly from 
higher modes of those loops. For all of the observable loops (that is, less 
than tens of light years across), formation occurred during radiation era 
(that is, when the universe was less than 10^ years old). As in Hogan [IB] 
we use scaling to normalize to the high frequency results from R.R. Caldwell 
and B. Allen [H]: = 10-^(G/i/10-9)^/VH7«/5)^/^ where p=l. At 

the high frequency limit of our curves this relation is used to normalize the 
parameterization for the number of loops created at a given time. 

7.2.3 Loop Lengths 

Loop sizes are approximated by the "one-scale" model in which the length 
of created loops scale with the size of the Hubble radius. The infinite strings 
create loops with a size that scales with the expansion, creating a distribution 
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with a characteristic size a/H. 

The expression for the average length of newly formed loops is {L{t)) = 
acH{t)~^, where a measures the length of the loop as a fraction of the Hubble 
radius. It is assumed that the newly formed, stable orbiting loops form a 
range of lengths at a given time, and this range peaks near a. For the 
numerical calculations, and to facilitate comparison with previous work, we 
use a delta function for the length. Thus only one length of loop is created 
at a given time and this is the average length defined above. Because of the 
averaging introduced by the expansion, this introduces only a small error 
unless the distribution of sizes is large, on the order of several orders of 
magnitude or more pSj, which would no longer be a one scale model. As 
a check, models of loop formation with several loop sizes are also analyzed, 
and it is found that the larger loops tend to dominate the background over 
the smaller. 

The loops start to decay as soon as they are created at time tc, and are 
described by the equation: 

L(t„ t) = acH{Q-' - 7G/i^^. (7.7) 

(This expression is not valid as L ^ but is an adequate approximation. As 
the length approaches zero Eqn. 17. 71 becomes less accurate, but in this limit 
the loops are small and contribute only to the high frequency red noise region 
of the gravitational wave background. This region is very insensitive to the 
exact nature of the radiation so Eqn. 17.71 is accurate within the tolerances of 
the calculations.) 

Loops that form at a time tc decay and eventually disappear at a time td- 
This occurs when L(tc, td) = and from Eqn. 17.71 

= acHit^y^ -^Gjj^-^^, (7.8a) 

td = tc + -k^- r.8b 

7G/i 

At some time t = t^ there are no loops created before time t^, so they no 
longer contribute to the background of gravitational waves. This is discussed 
in more detail below. 
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7.2.4 Loop Number Density 

For the number density of strings we parameterize the number created within 
a Hubble volume at a given time by Nt/a. The newly formed loops redshift 
with the expansion of the horizon by the inverse cube of the scale factor 
a(t)-3: 



n(t,t) = — - 
a \ 



H{t)Y f a{t) 
W) 



(7.9) 



n(t, t') is the number density at time t as seen by an observer at time t', and 
the redshift has been incorporated into the function. In the Appendix there 
is a description of the calculation of the scale factor. Analyzing the high 
frequency background equation given by C.J. Hogan |46j from R.R. Caldwell 
and B. Allen [TB] the curves are normalized by setting A^^ = 8.111. 

At an observation time t' the total number density of loops is found by 
summing over previous times t excluding loops that have disappeared, 

t' 

n,{t')=Y,n{t,t'), (7.10) 

t=te 

where tg is the time before which all the loops have decayed when making an 
observation at time t'. The earliest time in our sum, t^, is found by solving 
Eqn.O 

t' = te + 7^^, 7.11 

at a given t' . Thus t(.{t') tells us that at some time t', all of the loops formed 
before t^it') have decayed and are not part of the sum. For large loops and 
light strings trends toward earlier in the history of the universe for a given 
t', allowing more loops to contribute longer over cosmic history. For the final 
calculations, which are described in detail in the next section, t' and the 
number density are summed over the age of the universe. 

As the loops decay the energy lost becomes gravitational radiation, which 
persists in the universe and redshifts with expansion. This gravitational wave 
energy is detectable as a stochastic background from the large number of 
emitting loops. 
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7.3 Gravitational Radiation from a Network 
of Loops 



7.3.1 Frequency of Radiation 

The frequency of radiation is determined by the length of the loop, and has 
normal modes of oscillation denoted by n. The frequency of gravitational 
radiation from a single loop is given by, 

„ / N 2nc 

where tc is the creation time of the loop and t is the time of observation. There 
is also redshifting of the frequency with expansion which must be accounted 
for as this strongly affects the shape of the background curves. The strings 
radiate most strongly at the fundamental frequency so n = 1 for most of 
the calculations. This also facilitates computation times and introduces only 
small errors in the amplitude of the background, which is shown in Results. 
Ultimately the higher modes do not greatly alter the results at the resolution 
of the computations. 



7.3.2 Gravitational Wave Energy 

A network of loops at a given time t' radiates gravitational wave energy per 
volume at the rate: 

^-P^ = ^Gy.\n,{t'). (7.13) 
dt 

Recall that ns{t') is the number density n{t,t') summed over previous times 
t at a time of observation t'; given by Eqn. 17.101 

To find the total gravitational wave energy density at the present time 
the energy density of gravitational waves must be found at previous times 
t' and then summed over cosmic history. Since the energy density is given 
by Eqn. 17. 131 the total density is found by summing over ns{t') for all times 
t' while redshifting with the scale factor, a~^. Thus the total gravitational 
energy per volume radiated by a network of string loops at the current age 
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of the universe is given by: 



Pgwituniv) = iGfi^c — -^ns{t')dt' , (7.14a) 

Jto CLytuniv) j 

Here tuniv is the current age of the universe and to is the earhest time of loop 
formation. The internal sum is over t and the overall integral is over t' . In 
practice we can make to small to observe higher frequency, H loops, earlier 
during the radiation era. 

Eqs. 17.141 do not have explicit frequency dependence, which is needed to 
calculate the background gravitational wave energy density spectrum flgw{f) 
determined by Eqn. 17.21 This requires that the gravitational energy density 
be found as a function of the frequency: the number density as a function 
of time is converted to a function of frequency at the present time. This is 
discussed in detail in Appendix D. 

For most calculations it is assumed each length of loop radiates at only 
one frequency, / = 2c/ L. The number density, as a function of frequency and 
time t', is redshifted and summed from the initial formation time of loops to 
to the present age of the universe: 

p,M) = iGph -^P^^n{f^j^ t')dt'. (7.15) 

Note that the frequency is also redshifted. 

We have checked the effects of including higher mode radiation from each 
loop. The time rate of change of the gravitational energy density of loops at 
time t' is found by summing over all frequencies fj and weighting with the 
power coefficients Pj: 



dpgw{,f',t') 



df 



J2PjGp'cn{f'.,t'). (7.16) 



This is then integrated over cosmic time and redshifted to give the current 
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energy density in gravitational waves: 




(7.17) 



The power coefficients Pj are functions of the mode and have the form 
Pj oc [1061 lis] . This includes the behavior of cusps and kinks, which 

also contribute an f~^^^ tail to the power spectrum. Depending upon the 
percentage of power in the fundamental mode, the sum to infinity is not 
needed. As an example, if 50% of the power is in the fundamental mode only 
the first six modes are allowed. In general it is found that the higher modes 
do not significantly effect the values of the background radiation. Details are 
given in Results. 

From the gravitational radiation density we can find the density spectrum 
readily by taking the derivative with respect to the log of the frequency and 
dividing by the critical density: 

^,M) = (7-18a) 

A convenient measure is to take VLg^jh? to eliminate uncertainties in Ho = 
H{tuniv) = lOO/i km s~^Mpc~^ from the final results. (Note here h is the 
value of the Hubble parameter, not the strain on spacetime.) 



7.4 Results for the Stochastic Background 
7.4.1 General Results 

All of the computations show the expected spectrum for loops: old and de- 
cayed "H" loops leave a flat spectrum at high frequencies and matter era "P" 
loops contribute to a broad peak, with smooth transition region connecting 
them that includes comparable contributions from both. The frequency of 
the peak is strongly dependent upon the string tension, with lighter strings 
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leading to a higher frequency peak. This is a result of the lighter loops decay- 
ing at a lower rate, so smaller loops survive for a Hubble time. This behavior 
is exhibited in Fig. 17.11 

Fig. 7.2 shows the peaks of power output per unit volume per log fre- 
quency as a function of frequency at different times in cosmic history. The 
most recent times contribute to the broad peak, while later times contribute 
to the fiat red noise region. Although the power output per volume is larger 
earlier, the summation time is foreshortened so the contributions overlap to 
form the lower amplitude fiat region in the background. 

The differences in these calculations from previous publications can be 
accounted for by the use of larger loops and lighter strings as predicted in [IS] . 
The change in cosmology, i.e. the inclusion of a cosmological constant, is not 
a cause of significant variation. Since the cosmological constant becomes 
dominant in recent times the increased redshift and change in production 
rate of loops does not strongly affect the background. 

7.4.2 Varying Lengths 

The effect of changing the size of the loops is shown in Fig. 3. A larger 
a leads to an overall increase in the amplitude of but a decrease in 
the amplitude of the peak relative to the fiat portion of the spectrum. As 
expected there is little if no change in the frequency of the peak. In general 
the red noise portion of the spectrum scales with a^/^ and this is seen in 
the high frequency limits of the curves in Fig. 3 [16]. An increase in a 
leads to longer strings and an increase in the time that a given string can 
radiate. This leads to an overall increase in amplitude without an effect on 
the frequency dependence of the spectrum peak. 

In Fig. 17.41 a smaller a is plotted and shows the reduction in the back- 
ground. The general shape of the background remains unchanged, except 
for some small differences in the heavier strings. For heavier strings the am- 
plitude of the peak tends to increase relative to the red noise portion of the 
background spectrum, and the frequency of the peak shifts less with G\x. 
These differences become more pronounced for heavy and very small strings; 
and at this limit our results are found to match very closely those computed 
in [H]. 
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Figure 7.1: The gravitational wave energy density per log frequency from 
cosmic strings is plotted as a function of frequency for various values of G/i, 
with a = 0.1 and 7 = 50. Note that current millisecond pulsar limits have 
excluded string tensions Gfi > 10~^ and LISA is sensitive to string tensions 
~ 10~^^ using the broadband Sagnac technique, shown by the bars just below 
the main LISA sensitivity curve. The dotted line is the Galactic white dwarf 
binary background (GCWD-WD) from A.J. Farmer and E.S. Phinney, and 
G. Nelemans, et al., [STJ [76]. The dashed lines are the optimistic (top) 
and pessimistic (bottom) plots for the extra- Galactic white dwarf binary 
backgrounds (XGCWD-WD) [ST]. Note the GCWD-WD eliminates the low 
frequency Sagnac improvements. With the binary confusion limits included, 
the limit on detectability of Gfi is estimated to be > 10^^^. 
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Figure 7.2: The power per unit volume of gravitational wave energy per log 
frequency from cosmic strings is measured at different times in cosmic history 
as a function of frequency for G/i = 10~^^, a = 0.1, and 7 = 50. Note that 
early in the universe the values were much larger, but the time radiating was 
short. The low frequency region is from current loops and the high frequency 
region is from loops that have decayed. This plot shows the origin of the peak 
in the spectrum at each epoch; the peaks from earlier epochs have smeared 
together to give the flat high frequency tail observed today. 
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Figure 7.3: The gravitational wave energy density per log frequency from 
cosmic strings is shown with varying a for G/i = 10~^^ and 7 = 50. Here a 
is given values of 0.1, 10~^, and 10~^ from top to bottom. Note the overall 
decrease in magnitude, but a slight increase in the relative height of the 
peak as a decreases; while the frequency remains unchanged. The density 
spectrum scales as a^/^, and larger loops dominate the spectrum over small 
loops. 
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Figure 7.4: The gravitational wave energy density per log frequency from 
cosmic strings for a = 10~^ and 7 = 50. Smaller a reduces the background 
for equivalent string tensions, with a LISA sensitivity limit of Gn > 10~^^ 
if the Sagnac sensitivity limits are used. The bars below the main LISA 
sensitivity curve are the Sagnac improvements to the sensitivity, while the 
dotted line is the Galactic white dwarf binary background (GCWD-WD) 
from A.J. Farmer and E.S. Phinney, and G. Nelemans, et al., |37l[76]. The 
dashed lines are the optimistic (top) and pessimistic (bottom) plots for the 
extra-Galactic white dwarf binary backgrounds (XGCWD-WD) [27]. Note 
the GCWD-WD eliminates the low frequency Sagnac improvements, and 
increases the minimum detectable Gfi to > 10~^^. G/j, ranges from the top 
curve with a value 10"'' to the bottom curve with 10~^^. 
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7.4.3 Luminosity 

A residual factor of the order of unity arises from uncertainty in the typical 
radiation losses from loops when all modes are included averaged over an 
entire population. Changing the dimensionless luminosity 7 affects the curves 
as shown in Fig. I7.5[ Larger 7 decreases the amplitude of the red noise region 
of the spectrum, but increases the lower frequency region. Smaller 7 increases 
the amplitude of the peak and the flat portion of spectrum, while leading 
to a decrease in the high frequency region. This occurs because the loops 
decay more quickly for larger 7 and have less time to contribute over the 
long time periods of cosmology. In the matter dominated era the higher 
luminosity allows the loops to contribute more to the high frequency region, 
although it does not increase the peak amplitude. Higher luminosity also 
decreases the frequency of the overall peak by contributing more power in 
current loops. The power and change in length depend on the luminosity 
by -E" oc 7 and L oc —7. Note that the gravitational wave power depends 
on the square of the mass density but the first power of the luminosity, thus 
different behavior is expected varying 7 or /i. An analytic description of the 
scaling of the amplitude and frequency of Vtgw with 7 is given in [16] . 

7.4.4 Multiple Lengths 

Some models indicate several characteristic scales for the loop populations, 
with some combination of several scales. This implies multiple distributions 
of lengths, which can be described as a combination of different a's. Fig- 
ures 17. H 17.31 and 17.41 indicate the eventual dominance of the larger loops 
over the smaller in the final spectrum. Even though, in this model, the 
number of loops created goes as a~^, their contribution to the background 
is still much smaller and scales as a^/"^ . Thus an admixture of loop scales 
tends to be dominated by the largest stable radiating loops, which is shown 
in Fig. [721 

In Fig. I7.6l we assume 90% of the loops shatter into smaller sizes a = 10^^ 
while 10% are larger with a = 0.1. A "two scale" model is used and two 
relatively narrow distributions are assumed, as in the one scale model. It is 
imagined that the inital loops are large and splinter into smaller loops, with 
the percentage that shatter depedent upon the model. Overall the amplitude 
decreases by an order of magnitude below a pure a = 0.1 amplitude, but it 
is an order of magnitude above the a = 10~^ spectra. Thus the large loops, 
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Figure 7.5: The energy density in gravitational waves per log frequency from 
cosmic strings with varying luminosity 7 for G/x = 10~^^ and a = 0.1. Here 
the 7's are given as follows at high frequencies: 50 for top curve, 75 for 
the middle curve, and 100 for the bottom curve. The peak frequency, peak 
amplitude, and amplitude of flat spectrum all vary with the dimensionless 
gravitational wave luminosity 7. 
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Figure 7.6: The gravitational wave energy density per log frequency from 
cosmic strings for 7=50, Gfi = 10^^^, and a combination of two lengths 
(a's). 90% of the loops shatter into loops of size a = 10~^ and 10% remain 
at the larger size a = 0.1. The overall decrease is an order of magnitude 
compared to just large alpha alone. From bottom to top the values of G/i 
run from 10^^^ to 10^^. 
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even though they form a small fraction, contribute more to the background. 
Because of this weighting a one-scale model is likely to be suitable if appro- 
priately normalized. Note that a string simulation requires a dynamic range 
of scale greater than a to make a good estimate of a and to choose the right 
normalization. 

7.4.5 High Frequency Modes, Kinks, and Cusps 

Calculations are done which include higher mode behavior of the loops; this 
sends some of the gravitational wave energy into higher frequencies for each 
loop. Adding higher modes tends to smooth out, to a larger degree, the curves 
and slightly lower the peaks on the spectra. When applied to the power per 
volume per log frequency as in Fig. 17.21 we see a less steep high frequency tail 
at each time and a slight decrease in the amplitude of the peak. From Fig. 17.71 
a comparison of the fundamental mode versus the fundamental mode with 
higher modes included shows: a) the peak decreases, b) the peaked portion 
widens slightly, and c) the red noise amplitude remains relatively unchanged. 

Fig. 17.71 demonstrates a calculation in which ~ 50% of the power is in 
the fundamental mode and the power coefficients go as in the higher 

modes. For the fundamental mode Pi = 25.074, and the next five modes 
are = 9.95, P3 = 5.795, P4 = 3.95, P5 = 2.93, and Pg = 2.30. Their 
sum is ^ P„ = 7 = 50. The amplitude of the peak shifts from 1.202 x 
10~^ to 1.122 X 10~^, a 6.7% decrease. The frequency of the peak shifts 
more substantially, from 3.162 x 10~^ Hz to 7.080 x 10~^ Hz, an increase of 
124%. This increase in frequency of the peak does not affect the limits of 
detection by LISA because the amplitude shrinks very little. Note that the 
high frequency red noise region is essentially unchanged, and for Gfi = 10~^^ 
this is the region in which LISA is sensitive. 

Plots of spectra for other Gfi with high frequency modes show changes 
similar to those in Fig. 17. 7[ The peak amplitude drops slightly and the 
frequency of the peak is increased. Overall, the limits from Table [7l2] remain 
unchanged by the introduction of 50% of the power in the higher frequency 
modes. 

Cusps and kinks in the loops are responsible for high frequency "bursts" 
of gravitational wave energy [271 l2Hl [29l |88] . They are two different manifes- 
tations of string behavior: cusps are catastrophes that occur approximately 
once every oscillation, emitting a directed burst of gravitational radiation, 
while kinks are small wrinkles which propagate on the loops, emitting a 
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Figure 7.7: The gravitational wave energy density per log frequency for 7=50, 
Gfi = 10"^^, a = 0.1; plotted in one case (dashed) with only the fundamental 
mode, while for another case (solid) with the first six modes. The dashed 
curve has all power radiated in the fundamental mode, while the solid curve 
combines the fundamental with the next five modes: Pi = 25.074, P2 = 
9.95, P3 = 5.795, P4 = 3.95, P5 = 2.93, Pg = 2.30. Note the slight decrease 
in the peak and its shift to a higher frequency. The LISA sensitivity is added 
as a reference. 
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higher frequency of directed gravitational energy. The hght strings make it 
less likely to detect these random bursts; in particular they are more difficult 
to pick out from their own confusion background [l6]. Since cusps, kinks, 
and higher modes add an f~^^^ tail to the power spectrum of each loop, it is 
reasonable to expect that they will have a similar effect on Qgw The general 
behavior of high frequency contributions is represented in Fig. 17.71 

7.5 Millisecond Pulsar Limits 
7.5.1 Limits on Background 

Millisecond pulsars have been used as sources to indirectly measure a gravi- 
tational radiation background for a number of years [321 1201 ESI EZl [S21 172] . 
Limits on the strain, and correspondingly the energy density spectrum flgw, 
have been estimated recently at three frequencies: 1/20 yr~^, 1/8 yr~^, and 
1 yr~^ |53]. The limits on the strain vary with the frequency dependence 
assumed for the source of gravitational radiation at the corresponding fre- 
quencies: 

h,{f)=A(J-^ , (7.19) 

where P is the power of the frequency dependence of the source and is denoted 
as a in F.A. Jenet, et al.,[S3]- A table of values is provided in that source 
for various cosmological sources of stochastic gravitational radiation, with 
f3 = —7/6 used for cosmic strings. From [53], the cosmic string limits are: 

^gwil yr'^)h^ < 9.6 X 10"^ (7.20a) 
fig^(l/8 yr~^)h^ < 1.9 x 10~^ (7.20b) 
VLg^{l/20 yr^^)h^ < 2.6 x 10"^ (7.20c) 

Results of our calculations show that for heavier strings at the measured 
frequencies f3 is approximately -1, but for lighter strings it varies from this 
value. The variation of (3 is summarized in Table I7.H which has been calcu- 
lated at / = 1/20 yr"^ and / = 1 yr~^. This change in (3 will result in limits 
differing slightly only at low frequencies from those in Eqs. 17.201 [53] . 

In Figures 17.11 and 17.41 the limits for millisecond pulsars [53] are shown, 
with allowed values of Qg^ indicated by the direction of the arrow. For large 
strings a = 0.1, the maximum allowed Gfi is about 10~^, Fig. l7.1[ Decreasing 
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Table 7.1: /? for given G/i at / = l/20?/r ^ and f = 1 yr ^, a = 0.1, 7 = 50, 
wliere he (x f^. 







P 




f = 1/20 yr-i 


f = 1 yr'^ 


10-6 


-1.00306 


-1.00031 


10-8 


-1.05314 


-1.00894 


10-9 


-1.05985 


-1.03853 


10-10 


-0.931749 


-1.06593 


10-11 


-0.735778 


-0.987882 


10-12 


-0.682874 


-0.778679 


10-13 


-0.676666 


-0.688228 


10-14 


-0.676034 


-0.676054 


10-15 


-0.675971 


-0.674794 


10-16 


-0.675965 


-0.674668 



a decreases tlie amplitude of the background allowing for heavier strings 
within the pulsar measurement limits. A more complete list is displayed in 
Table El 

7.5.2 Dependence of Characteristic Strain on Frequency 

From the curves of Fig. 17.11 one can find the frequency dependence of the 
characteristic strain as given by hc{f) = A{f /yr~^)^ . Again, (3 is the power 
law dependence of he on frequency, and is generally written as a with a 
numerical value of -7/6 [SSj EHl EH]- In this paper /? is used to prevent 
confusion with the length of loops. If Qg^ is known, then using Eqn. 17.23) we 
find: 

27^2^2 pil+13) 

^gM) = 3ij2 (3.17 X 10-8 ('^•^^^ 

With this equality P is found by computing dhi{flgyj)/d\ia f at / = lyr^^ = 
3.17 X lO^^Hz and setting this equal to 2(1 + /3). The results are listed in 
Table O 

From Fig. 17.11 the change in slope with Gfi is evident at the measured 
frequency fo = 1 yr~^. For large values of Gfj,, fo is in the fiat red noise 
section of the spectrum so the slope of flgw goes to zero, and /3 approaches 
-1. As the string tension is reduced, /3 is measured first at the high frequency 
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tail so the slope is negative, and thus (3 is more negative. Once over the 
peak on the high frequency side, the slope becomes positive, thus /? must 
be greater than -1. The characteristic strain power law dependence j3 then 
approaches a value of -0.67. 

This dependence on G^x changes the value of the constant A in Eqn. 17.191 
and thus the limits on VLg^ [33]. Higher mode dependence only modestly 
changes the values of (3 since it approaches the limits -1 and -2/3 for high 
and low frequencies respectively. The frequency at which it takes on these 
values depends upon G/i as seen in Fig. 17.11 



7.6 Future LISA Sensitivity 
7.6.1 Physical Principles 

The pending launch of LISA will give new opportunities to observe gravita- 
tional waves in general, and a stochastic background in particular [221 HI 1231 
[To] . There are a number of theoretical sources for this background, including 
cosmic strings [Ml EH EHl ESI HT] . LISA sensitivity has been modeled and 
methods for improving sensitivity by calibrating noise and integrating over a 
broad band have been proposed [621 HHl [2l] . The cosmic string background 
density spectrum calculated here is compared to the potential for detection 
by LISA. 

Detectors such as LISA measure the strain on spacetime, so a relation 
between the background spectrum and strain is needed. The rms strain is 
denoted by Kmsif) or h{f), and for an isotropic, unpolarized, stationary 
background [531 SHI [Ml H] ■ This causes problems since typically hrms is the 
average value of hij{x) over several wavelengths divided by 2. It is better to 
use the spectral density Sh{f), see Appendix D for details of these relation- 
ships. 

47r2 

^,M) = ^fSnif). (7.22) 

Another measure is the characteristic strain spectrum hc{fY = 2fSh{f), 
when substituted gives: 

^.4/) = ^J'hl{f). (7.23) 
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In the literature, Sh{f) is often used and referred to as the spectral density 
[69]. These relations are important as they relate the energy density spectrum 
^gw{f), with the strain that LISA will detect. This limit on the detectable 
strain is shown as a curve on our plots and is taken from j62j . 

Here the LISA sensitivity curve is calculated with the following parame- 
ters: SNR=1.0, arm length=5xl0^ m, optics diameter=0.3 m, laser wavelength=1064 
nm, laser power=1.0 W, optical train efficiency=0.3, acceleration noise=3 x 
ni/(s2\/ffi), position noise=2 x 10~^^ m/y/Hz. 

Extracting the background signal by combining signals over a broad band 
and an extended time can increase the sensitivity of the detector. One 
method known as the "Sagnac technique" is shown on plots of the back- 
ground [ISl Ell, for example see Fig. 17.11 

7.6.2 Limits for LISA Sensitivity 

For both Figures 17.11 and 17.41 there are limits on the detectability of the 
background due to strings of various mass densities by LISA from [52] • From 
Fig. 17.11 with a = 0.1 the minimum value of the string tension measurable 
is ~ 10~^^. In Fig. 7.4, which plots a smaller a = 10~^, the minimum 
is G/i ~ 10^^^. These values and others are listed in Table 17.21 Sagnac 
observables increase our available spectrum and change the detectable string 
tensions as indicated ^Hj . 

Shown in Table 17.21 are the limits on string tensions for string sizes a 
of 0.1 to 10~^, along with the corresponding minimum tension that LISA 
can detect. Galactic and extra-Galactic binaries confusion limit is taken 
into account for this table. For a < 10^^ the maximum string tension will 
be greater than 10~^, although these values are excluded by WMAP and 
SDSS |109j . It is clear from the table and figures that decreasing a increases 
the maximum string tension allowed as well as increasing the minimum string 
tension detectable by LISA. 

7.6.3 Confusion Noise from Binary Systems 

Galactic white dwarf binaries add a background "confusion" that shrouds 
portions of LISA's spectrum. This background has been studied and predic- 
tions given for its spectrum [76]. There are also extra- Galactic binary sources 
which similarly create confusion noise [57]. These confusion limits are shown 
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Table 7.2: Millisecond Pulsar Limits and LISA sensitivity for 7 = 50. WMAP 



and SPSS [T09j have excluded G/i > 10'^. 





Millisecond Pulsar 


Minimum LISA 


a. 


Limits, G/i < 


Sensitivity, G/i > 


0.1 


10-9 


10-16 


10-2 


10-8 


10-15 


10-3 




10-14 


10-^ 




10-13 


10-5 




10-12 


10-6 




10-11 



on plots of the gravitational radiation spectra, and they affect the minimum 
values of string tension that are detectable. 

7.7 Conclusions 

Broadly speaking, the results here confirm with greater precision and reli- 
ability the estimates given in [36]. Our survey of spectra provides a set of 
observational targets for current and future projects able to detect stochastic 
gravitational wave backgrounds. For light strings, the string tension signif- 
icantly affects the spectrum of gravitational radiation from cosmic string 
loops. In addition to a reduction in the overall radiation density, the most 
significant change is the shift of the peak of the spectra to higher frequen- 
cies for light strings. Near the current pulsar limits, the peak happens to lie 
close to the frequency range probed by the pulsars, so spectra as computed 
here are needed as an input prior for establishing the value of the limits 
themselves. For much lighter strings, the peak is at too high a frequency for 
strings to be detectable by pulsar timing as the spectrum falls below other 
confusion backgrounds. Lighter strings will be detectable by LISA. For the 
lightest detectable strings the peak happens to lie in the LISA band and the 
detailed spectrum must again be included in mounting observational tests. 
The spectrum from strings is quite distinctive and not similar to any other 
known source. 

A high probability of forming loops with stable, non-self-intersecting or- 
bits, as suggested by recent simulations, leads to larger string loops at for- 
mation which in turn give an increased output of gravitational wave energy 
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and improved possibilities for detection for very light strings. Recent sim- 
ulations have shown stable radiating loops form at a significant fraction of 
the horizon, of order a = 0.1; for loops of this size our calculations show 
that the maximum string tension allowed by current millisecond pulsar mea- 
surements is G/j, < 10~^, and the minimum value detectable by LISA above 
estimated confusion noise is G/x ~ 10~^^. In field theories, the string tension 
is typically related to the scale by Gfi oc A^/m^^, so the maximum detectable 
scale in Planck masses currently allowed by millisecond pulsars is A < 10^^'^, 
or around the Grand Unification scale; with LISA, the limit will be about 
As ~ 10~^ or 10^^ GeV. These results suggest that gravitational wave back- 
grounds are already a uniquely deep probe into the phenomenology of Grand 
Unification and string theory cosmology, and will become more powerful in 
the future. The most important step in establishing these arguments for an 
unambiguous limit on fundamental theories will be a reliable quantitative 
calculation of a. 
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Chapter 8 
Conclusion 



The results of this study indicate that cosmic strings may have a strong 
gravitational wave signal from loops in the millisecond frequency range, cor- 
responding to the heart of the detection band for the Laser Interferome- 
ter Space Anetenna (LISA). This signal comes from both individual cosmic 
string loops within the Milky Way and a stochastic background from past 
and present loops throughout the universe. For light strings the galactic 
loops produce a large average signal that acts as a "confusion noise" . 

Observations from WMAP, Sloan Digital Sky Survey (SDSS), and mil- 
lisecond pulsars have precluded heavier strings, so this study has extended 
the range of previous investigations and undertaken some new ones. We 
have extended the range over previous investigations into the stochastic back- 
ground of gravitaitonal waves, and have done a new investigation into the 
detectability of the gravitaitonal radiation from the normal mode oscillations 
of cosmic string loops. 

From our calculations of the stochastic background it is found that for 
large loops string tensions down to G^i ~ 10~^^ may be detectable, see 
Fig. 17. 1[ For these large loops the upper limit is also set at G/i < 10^^ 
by analysis of millisecond pulsars. This further constrains string formation 
scenarios below the Planck scale. In field theory the energy scale of the phase 
transition is given by = Gixrn^i where mpi = 1.22 x 10^^ GeV is the Planck 
mass. From our large loops we have and upper bound of A^ < 10~^'^mpi and 
LISA has a sensitivity down to A^ ^ lO^^nipi. 

For the loops within the galaxy, our calculations indicate that light and 
large strings may have a significant signal in bulk, see Fig 16.31 The signal 
from this confusion noise from large loops is detectable by LISA down to a 
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string tension of Gjj, = 10~^^, which is lower than previous esimates. The 
energy scale detectable is down to = 10~^'^mp/ or 10^'^ GeV. 

Thus cosmic strings can probe fundamental theories of physics from en- 
ergies of 10^^-5 GeV to lO^-^ GeV {IQ-^-^nipi to IQ-^-^nipi). Theories with 
energy scales in this range include: 

• GUT scale: 10^^ GeV=lQ-^mpi 

• Peccei-Quinn scale: 10^ GeV = 10~^^mpi 

• B-L symmetry breaking scale: some theories 10^ to 10^° GeV (10~^^ to 
lO-^rUpi) 

• Superstring brane cosmologies: 10~^ < G/i < 10^^^ 

Thus the gravitational wave signal from cosmic string loops is a unique 
probe of field theory and string theory. Further refinement in the value of 
the loop size a is required to eliminate uncertainties in these results. 
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Appendix A 
General Relativity 



In this appendix we write out some of the basic results from general relativity. 
A more complete description can be found in any general relativity text [81 

[TnnimEn] 

First the Einstein-Hilbert action: 

Se = -J— I R (fx (A.l) 



where R is the Ricci scalar and g is the trace of the metric. We define the 
Christoffel symbol by, 

rj. = ^^7^" {d,g^u + d,g^^ - d^g^,) (A.2) 
and the Riemann tensor, 

The Ricci tensor is R^^i, = R'^^^, and the Ricci scalar is given by i? = g^^R^y. 
We define the matter action as: 

J4 



Sm= T^^'-gpud'x, (A.4) 



where T^'^ is the stress-energy (or momentum-energy) tensor. The stress- 
energy tensor measures the energy, energy flux, momentum density, and 
stress at a point in spacetime. The components are given by: Too= en- 
ergy density, Toj=energy flux, Tjo= momentum density, Tjj= pressure, and 
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Tij= shear stress. In matrix form: 



T, 



( ^ 


energy flux 




P. 


shear stress 


p 




Py 


\ density 


shear stress 


Pz 



(A.5) 



Thus the total gravitational action is the sum of the matter and Einstein 
actions, Sg = Se + Sm- Varying this action with respect to the space-time 
metric, gfj,i, and setting it to zero, 6S = we find: 



R 



which are the Einstein equations. 

If we take the trace of Eq. IA.6I we get R 
Eq. IA.6I gives us, 

„ „.„ 1 



(A.6) 

StiGT. Reinserting into 
(A.7) 



Eqs. IA.6I and IA.7I are the same, but Eq. IA.7I is more useful if the stress- 
energy tensor is known, as in the analysis of cosmic expansion. 
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Appendix B 
Cosmic Expansion 



In this appendix we review some basic results pertaining to the expansion of 
the universe, perhaps most importantly the dependence of the scale factor, 
a{t), with cosmic time. A more complete discussion can be found in any 
good text on the subject of cosmology [33l[60l[78l|73l[71l[ini[65l[79]. 



B.l Homogeneous and Isotropic Universe 

The space-time interval is defined by 

ds^ = g^^dx^dx", (B.l) 

where g^u is the metric of the space-time. 

For a homogeneous and isotropic universe the Friedmann-Roberston- Walker 
(FRW) is a suitable metric which defines the spacetime interval, and is given 
by: 

ds^ = -de + a{tf ( - r'dnA , (B.2) 
\1 — kr J 

where t is the cosmic time, r is a radial coordinate distance, k is the curvature, 
and a{t) is the scale factor. 

The coordinate distance vector f{t) is given by, 

f{t) = a(t)e, (B.3) 

— * 

where ^ is a comoving coordinate that remains constant with the expansion. 
Thus a{t) defines the relative size and rate of expansion of the universe, and 
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can be solved using the Einstein equations if the constituents of the universe 
are known, and thus the stress-energy tensor is known. 

The Hubble law states that the velocity of a galaxy is proportional to its 
distance from the observer: 

f = - r, (B.4) 
a 

= Hr, (B.5) 

where H is the Hubble parameter and varies with cosmic time. The Hubble 
law is not valid for very large redshifts, where cosmic expansion changes due 
to a transition in the dominant constituent of mass/energy. 

Now we define the constituents of the universe and determine their effect 
on the expansion, and therefore calculate the scale factor, a{t). Einstein's 
equations, 

Rf.u = 8nG (t,,, - ^g^,T^ . (B.6) 

are used with the FRW metric. 

Lets start with an isotropic ideal fluid, with the stress-energy given by, 

T^u = (-P + p)u^Uu + Pg^,u, (B.7) 

where the are the four velocity of the fluid with respect to the observer. If 
our reference frame is at rest with respect to the fluid, that is = (1, 0, 0, 0) 
and locally flat we can write: 

T^. = iP + p)6o^6ou + Pv^u (B.8) 

where 77^1^ = (— 1,1,1, 1) is the Minkowski metric. Taking the trace to use in 
Eq.lXTn 

tr(T^,) = r = T;i = r7'^%„ (B.9) 
= -P + 3P. (B.IO) 

B.2 The Friedmann Equations 

After a bit of calculating one finds the Friedmann equations: 

a AtiG , ^ s 

- = -^3P + p, B.12 
a 6 
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where p is the sum of all mass/energy densities. These equations describe the 
dynamics of an isotropic and homogeneous universe, with the first equation 
usually referred to as THE Friedmann equation. 

This first equation is often rewritten in terms of the "criticial" density, 
the density of the universe which results in a curvature of zero (A;=0). This 
is given by. 

Where Hq is the current value of the Hubble parameter. We further define 
the densities in terms of the critical density, 

= ^, (B.14) 

Pc 

where the j represent the various components of the cosmic soup. This leads 
to a rewrite of the Friedmann equation, 

- HlY,^, - (B.15) 

3 

Observations from WMAP have placed ilneto — 1.02 ±.02 and thus A; 0. 
Regardless, the density of the universe is startlingly close to unity, a fact 
which still requires much explanation to be understood. 



B.3 Cosmic Constituents 

Now we introduce a small bestiary of the cosmos. These components gener- 
ally have densities that vary with temperature, and thus with the expansion. 
In general we can write down the equation of state of each component in 
terms of the energy density, 

P = wp (B.16) 

where w is a scalar. The reader should be careful to discern the difference 
between the momentum and pressure. As is the convention the same variable 
p is used for both. For this chapter we will use the upper case for the pressure. 

Conservation of the stress-energy tensor gives a relationship between the 
scale factor and density, which in the curved FRW spacetime must be written 
in covariant form: 

D^T^"" = 0, (B.17) 
d^T^''' + r'^^T^'' + r''^^T^'^ = 0, (B.18) 
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where F^^ are the Christoffel symbols. This can also be written in the 
form: 

-^9^v^T'^^ = 0, (B.19) 

where g is the determinant of the spacetime metric. 

Using the = component of the stress-energy tensor leads to, 

p + 3{p + P)- = 0. (B.20) 

a 

This equation gives the time dependance of the density of each constituent, 
and also describes the dependance of the densities on the scale factor (or 
volume). 

The constituents of the universe include radiation (utra-relativistic), mat- 
ter (non-relativistic), vacuum or "dark" energy, neutrinos, gravitational waves, 
and cosmic strings of the infinite and loop variety. It has been found via 
observation that gravitational waves and cosmic strings form a miniscule 
fraction of the cosmic energy density, and thus aren't considered when cal- 
culating the scale factor. Since neutrinos have mass, it is likely that some of 
the species are non-relativistic at the present time. This is a minimal con- 
sideration given that their contribution is at least five orders of magnitude 
smaller than the dark energy and matter at the present. See Eqs. ID. 91 for 
values of the density parameters, fl of the various components. 

To find the equations of state of the components we start with the en- 
ergy density by integrating over all possible microstates assuming a non- 
interacting fiuid: 

^ f 9 d^P e 

Aug Ve^ - de 

(2^L e^±l 

where = + is the energy of each particle, g is the number of internal 
degrees of freedom, and the plus sign is for fermions and minus for bosons. 
For a non-relativistic gas the the energy density is approximately just given 
by the mass density times c^. For an ultra-relativistc gas m = and the 
results simplifies greatly. 



(B.21) 

(B.22) 
(B.23) 
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To find the pressure we must sum over all particles in a hemisphere that 
collide with some area AA n. These particles then impart mometum to the 
surface given by 2(p ■ n) for a total contribution to the pressure given by 

2(p-n)AN 

AP = ' , B.24 

AtAA ^ ' 



( — de 



27re V ± 1 



3e \ e kT :izl 



cos2(^) sm{e)ded(j), (B.25) 

(B.26) 



The expression in parenthesis is the number density of particles with energy 
e to e + de. 

To find the total pressure we integrate over the energy states, 



2 TOO 2 ^2 

de, (B.27) 



m 3e e kT ±1 



P "^'^ -^t^" de. (B.28) 

Jm e fcT ±1 



3 67r2 



Immediately we see that massless (ultra-relativistic) particles have an equa- 
tion of state given by P = p/3 regardless of the particle species or chemical 
potential. 

For the equations of state for the constituents of the universe we have: 

Non-relativistic matter, ti? = 0, p oc a~^, (B.29) 
Ultra-relativistic matter, = p oc (B.30) 

Vacuum energy (Dark energy), w = —1, p = constant. (B.31) 

The results for the scale factor were found from Eq. IB.20[ These relations 
indicate a universe composed of "dark enegy" will eventually have the expan- 
sion dominated by it. Both relativistic and non-relativistic matter density 
shrink as the universe expands, which corresponds to the observation that 
the very early universe was dominated by radiation, which was replaced by 
non-relativistic matter, finally to be overtaken today by vacuum energy. 

A description of calculation of the scale factor for this work is found in 
Appendix D. 
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B.4 Entropy of the Universe 



We also define the entropy by, 

S^sa{tf, (B.32) 

where s is the entropy density. From the standard cosmological model ho- 
mogeneity implies that entropy is conserved as there are no temperature 
differences to drive heat flow, thus dS/dt — 0. 
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Appendix C 
Signal Detection 



In this appendix we characterize the signal on a detector and the mathe- 
matical fomalism used in this process. From an individual detector we can 

characterize the output by s{t) = h{t) + n{t) where h{t) = F^h+ + Fxhx 
is the gravitational wave signal in terms of its two polarizations, and n{t) is 
intrinsic noise. In order to extract the noise and to discern the signal in both 
space and time, we take the following steps. 



C.l Spectral Density and Characteristic Strain 

In the transverse traceless gauge we can define the wave in terms of its fourier 
components, in TT gauge by, 

/oo r p2n 

df dk dip {h+{f,k,ip)elj{k,2p) + 
■oo Jci Jo 

+ hM,ki^)e^ik,i^))e-''''^^'-'-''\ (C.l) 

where, hA{—f,k,'4') = h\{f,k,ip). Also ^ is a unit vector pointed in the 
direction of propagation and ■0 is the angle of polarization of the wave. The 
two polarization tensors have the orthogonality relation ef-(k)e^,(k) — 25^^ , 
no sum over A. Now we define the spectral density Sh{f) from the ensemble 
average of the fourier amphtudes, 

h*M, k i^yhA'if, k', i^')) = 5AA'5{f - f) 
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where S\k, k') = (5(0 - (j)')S{cos{e) - cos(0'))- 

Now we write down the average of the real waves we get, 



2 

rms ' 
oo 



dfS.if), 



f=oo 



d{\nf)fS,{f). 



(C.3) 
(C.4) 

(C.5) 



The characteristic strain is defined by, 



(/i,,(t,x)/i^^(t,x)> = 2 / dilnf)hl{f). 



/=oo 



/=0 



Therefore, 



hl{f) = 2fS,{f). 



(C.6) 



(C.7) 



From above we see that Sh{f) gives the rms value of h within a range of / if 
the strain is confined to a small region of /: 



h' ^ Af S,if), 
or more broadly for the Fourier transform. 



h{f) ^TS,{f). 



where T is the total integration time of observation. 
For LISA the Fourier transform is defined as [631. 



(C.9) 



Hf) 



which results in the equality. 



T/2 



VT J -t/2 



dt h{t)e'^^^\ 



(C.IO) 



Skif) = Hf) 



(C.ll) 



In this relationship, the spectral density is often defined by the new variable, 
Skif) = h). 
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C.2 Noise and Strain Sensitivity 

The fourier components of the noise satisfy, 

inKmAf')) = l^if - n6.,s^:\f), (c.12) 

where S'„(/) is the spectral noise density, dimensions Hz^^ So we get, 

POO 

(n'it)) = / dfSnif). (C.13) 
Jo 

We can also define the strain sensitivity as a/ Sn{f)- 

C.3 Single Detector Response 

The response from detetor is given by S{t) = h[t) + n{t) where h is the signal 
from the waves and n is the noise. 

h{t) = D''h,,{t), (C.14) 

where D^^ is the detector tensor. Examples include an interferometer given 
by D^^ = ^{rhiTTij — nifij), and a bar detector Dij = iiij. Note that m and n 
are unit vectors along the axes of the arms, not necessarily orthogonal. From 
eqn. IC.ll we get, 

/oo p /»2-7r 

df dk # 

X hA{f,k,^P)D'^et^{k,^)e-''^^^'~'■'\ (C.15) 
We can also define the detector pattern functions as, 

FACk,i^)=D''et^{k,i^), (C.16) 

therefore, 

/CO p p27r 

df dk di) 
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and its Fourier transform, 



X h4f, k, ip)FA{k, ,py^-f^-\ (C.18) 

For a stochastic background we find that < h{t) >= so we want to look 
at < h?{t) >, 

POD 

{h'{t)) = F / dfSnif), (C.19) 
^ " Itr^^ f^AV.)F-AW. (C.20) 

j4 — -|-, X 

where F is the pattern function of the detector and Sh{f) is the spectral 
density. These equations are ultimately independent of ip for a stochastic 
background, so when calculating backgrounds they are omitted in future 
formulae. 

C.4 Signal to Noise 

A gravitational wave of Fourier transform h{f) incident on a dector has a 
signal to noise ratio of, 

(57Vi?)2 = 4j^"rf/M^. (C.21) 
This is an upper limit assuming optimal filtering. 

C.4.1 Signal to Noise of a Single Detector for a Peri- 
odic Source 

A period source delivers a signal at the detector given by, 

h{t) = F+{k,i^)h+{t) + F^{k,i^)hAt), (C.22) 

where hA{t) are periodic functions of time with some frequency /„. The two 
/i's may differ up to a phase (p. 
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Prom the equation above we find, 



oo 



F+{k,i^)h+,o + F^{k,i^)h^,oe"^' (C.24) 
where T is the total observation time. 



C.4.2 Signal to Noise of a Single Detector for the Stochas- 
tic Background 

To compute the signal to noise ratio we need both {h'^it)) and {n^{t)). From 
equations above we note that: 

{s\t)) = {n\t)) {h\t)) , (C.25) 

/•oo 

= / df{SM + FSn{f)). (C.26) 

To calculate the signal to noise the signals must be put into bins corre- 
sponding to the observation time T of frequency range A/ = 1/T. We then 
sum over the bins such that J Sndf = ^^^A/, and look in each freqency 
bin to find the signal to noise (SNR): 

SNR = (C.27) 

- (C.28) 



SnifW 
Snif) ■ 



(C.29) 



The final result is independent of integration time. Thus for a single de- 
tector searching for a stochastic background signal, it will either stand out 
immediately or cannot be detected over the noise. 

The minimum detectable value of Sh{f) for a given SNR and Sn{f) is 
given by, 

{SniDUn = ^SNR. (C.30) 
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The corresponding minimim detectable spectral density D,guj is then found to 
be, 

{^,M)Un = ^^SNR. (C.31) 

C.5 Two-detector Correlations 

Now we want to look into multiple detectors, and calculate the correlations 
between these detectors. For the i^^ detector we get the signal, 

/oo r 
df / dk JiAif, k)F^^ik)e-''''^^'-'-'^\ (C.32) 

SO we can write the fourier components of h{t), 

%)(/)= E [ dkhA{f,k)F^){k)e''-f'-^^% (C.33) 

So correlating two amplitudes from two different detectors will be given by, 

fT/2 

-r/2 

Q will tend to fall off quickly for large t — t', and the limit for large T, 



.T/2 . 

^12= / dt dt'si{t)s2{t')Q{t,t'). (C.34) 

J-T/2 J 



/oo 
dfsl{f)S2{f)Q{f). (C.35) 
-oo 

We want to find the ensemble average < 5*12 >, and the contribution from 
GW's is given by, 

{hi2)) = r df(hi;h^2))Q{f), (C.36) 

J —oo 

/oo p 
df / dkdk' e^'^'f ik-Si-k'-S2) J2 X 
-oo J 

F^{k)Fi{k') (h*M, k)hA'{f, k')) Qif), 
= T df^Sn{f)T{f)Q{f), (C.37) 
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where S{f — f) — dt — T and 



r(/) 



dk 
An 



■Ax 



(C.38) 



(C.39) 



C.5.1 Filtering 

We define the variation of 5*12 as, 

N = Sn- {Su) , 
which leaves < N >= 0. 

(N^) = {S!,)-{S,,)\ (C.40) 

dfdf'Q{f)Q%f') {si{f)S2{m{f')s;{f) 

Sl(m(f))(s;(f)Sr{f))]. (C.41) 



We are primarily interested in the regime where n{t) » h(t), and the 
noise in the two detectors is uncorrelated. Then Si ^ rii, and we also assume 
a gaussian distribution, 



rp poo 2 

{N')^- df Qif) Slif), 

^ J —oo 



where. 



Sn{f) = 4^^\f)s¥u)- 

For the signal to noise ratio we define. 



SNR = { 



1/2 



(C.42) 



(C.43) 



(C.44) 



Allen defines it without the square root. 

We are looking for a filter function Q{f) that maximizes the SNR. We 
can define an inner product using the conventions above as. 



dfA\f)B{f)Sl{f). 
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(C.45) 



Which gives us, 

(7V2> = |(g,g), (C.46) 

and 

(V2)> = |(4^) (C.47) 
Then we must maximize, 

,C.48) 



The solution is given by, 

r 

n(f \ = h— 

Slif) ' 



^ k^^^y^, (C.49) 



with k an arbitrary constant. 
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Appendix D 

Numerical Calculations for the 
Stochastic Background 

D.l General 

For computations the relevant differential and integral equations are solved 
numerically. Here is a list of the equations that are used, most of which have 
been discussed in the chapters above. In addition there are several more that 
are expanded upon in this appendix. Below note that n{t, t') is the number 
density with time t at an observation time t', while n{f'. f') is the number 
density as a function of frequency and time t' . The time t is always earlier 
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than the time t', which is the time of observation of the quantity in question. 



P = Prad + Pu + Pmatter + PA 

+ Poo + PL + P3«), (D-2) 

L{t,,t) = acHit,)-^ -GiJ—^, (D.3) 



PgwU) = iGfX^C j 

Jtr 



to Mn.™ 

Recent observations allow us to neglect loops pl, infinite strings poo, and 
gravitational waves pgw in the Friedmann equation, so these terms are re- 
moved [llOj . 



D.2 Scale Factor 

Calculations of the scale factor are done using the standard cosmological 
model with a homogeneous isotropic spacetime defined by the Friedmann- 
Robertson- Walker metric, solved using Einstein's equations [HlESlEn]- As an 
addition the effects of a cosmological constant on the cosmic string gravita- 
tional wave background are computed. The energy densities used to describe 
the dynamics of the scale factor are radiation, neutrinos, matter, and a cos- 
mological constant. The equation of motion for the scale factor is determined 
by the Friedmann equation. 

2 

d \ SvrG , \ A „N 

{Prel+Pm) + —, (D-7) 



3\r I ' r til' / ' o 
6 

where the p's are functions of the scale factor. We then integrate the above 
equation to find a{t). 
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The relativistic and matter densities follow the relations: 



P,„ ^ (D.8b) 

The gravitational radiation density also goes as a{t)~^. 

In order to compute the scale factor the Friedmann equation is trans- 
formed into an integral equation and solved numerically. The main source of 
error is the relativistic treatment of neutrinos for all time, but this will incur 
a very small error in more recent times. Given that neutrinos have mass, it is 
possible several species will be non-relativistic at the present time. This is a 
minimal consideration as the matter and dark energy contributions are larger 
by at least five orders of magnitude when the neutrinos are non-relativistic. 
The following parameters are taken from the Particle Data Group, 2006 [llOj: 



(D.9a) 

^matter = 0.24, (D.9b) 

(D.9c) 
(D.Qd) 
(D.9e) 

(D.9f) 

where Qj = Pj/pc, and the constant h is given by Ho = lOO/i km s-^Mpc-^. 





= 0.76, 


"matter 


= 0.24, 


^rad 


= 4.6 X 10"^ 




= 3.15 X 10"^ 


h 


= 0.73, 






pc 


SttG" 



D.3 Gravitational Radiation from Loops 

At a given time t', there is an earlier time t before which all of the loops have 
decayed, which we denote te- To find this value at any time t', Eqn. 17.71 is 
set to zero giving L{t, t') = 0, and t = is calculated numerically. Note that 
te is a function of t'. For very early times te{t') can be lower than the initial 
time in our sum over cosmic history to. In this case to was used instead of 
the value for t^it'). 

For many of the calculations it is assumed each loop radiates at only one 
frequency, / = 2c/ L. Then a numerical list of values is computed for t = t^ 
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to t — t', {L{t , t') , n{t , t')} , and an interpolation function is created from the 
list. This is the number density as a function of length and time t', n[L,t'). 
Next a list of frequency and number density is created, {2c/L, t')}, with 
lengths running from L{te,t') to L{t',t'). An interpolation function is made 
from this list, and this results in the number density as a function of time t', 
n{f',t'). But the frequencies must be redshifted to the current time tuniv so 
the function is n{fa{t')/a{tuniv)-,t'). 

Then we must integrate over all t' from the earliest time of loop formation 
to to the present age of the universe tuniv This integral, 

is transformed into a sum, 

pUf) = iGp'c ^(/#\' At'. (D.ll) 

Suitable At"s are chosen to ensure both computability and to minimize er- 
rors. In general we solve all equations in logarithmic bins so At' will vary 
in accordance. Using small values of At to facilitate computation speed will 
result in relatively small error; choosing At = 0.1 in log scale gives a 10% 
smaller value for Qgyj compared to At = 0.001 or At = 0.01 in log scale. 
Values smaller than At = 0.01 give little or no change to the background. 
After the sum the result is the gravitational wave density as a function of 
frequency at the present time. 

To include higher modes the power coefficients Pj must be included in the 
calculation. The subscript j is used instead of n here to prevent confusion 
with the number density, n{f',t'). At each time t' the number density, and 
thus the time rate of change of the energy density, of loops is calculated 
exactly as above, but more frequencies are included. Each number density 
at frequency fj is summed and weighted with Pj, and this number density is 
then integrated over cosmic history, 

J to (^yj'univ) 
1^ ^ij^univ) J 



(D.12) 
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The overall integral is then transformed into a sum as before. The sum of 
the power coefficients is = 7, where the Pj oc j^'^^'^. 

The final result is to calculate at the present time. This is done 

by finding the current energy density of gravitational waves as a function of 
frequency, and then taking the derivative with respect to f: 

U.Jf) = l^. (D.13) 

This is then plotted as \og{Qgyj{f) h"^) vs. \og{f/Hz), as seen in Figures 7.1, 
7.3 — 7.7. This quantity is, again, the density of stochastic gravitational wave 
energy per log of frequency in units of the critical density and is related to 
the strain on space by the gravity waves. 
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